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Abstract 

We formulate the laws governing the dynamics of a crystalline solid in 
which a continuous distribution of dislocations is present. Our formula- 
tion is based on new differential geometric concepts, which in particular 
relate to Lie groups. We then consider the static case, which describes 
crystalline bodies in equilibrium in free space. The mathematical problem 
in this case is the free minimization of an energy integral, and the asso- 
ciated Euler-Lagrangc equations constitute a nonlinear elliptic system of 
partial differential equations. We solve the problem in the simplest cases 
of interest. 
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Introduction and Summary 

Classical elasticity theory rests on the hypothesis that a solid may exist 
in a stress-free state in Euclidean space. The dynamics of the solid is then 
described by a time-dependent mapping taking the position of each material 
particle in this relaxed reference state to its actual position at a given time. 
Now the hypothesis of the existence of a relaxed state is violated when 
dislocations are present in the solid. These cause internal stresses even in 
the absence of external forces. 

Solids in which a continuous distribution of dislocations is present have 
been treated in the linear approximation ( |10j . [llj . |13j ). However, the 
conceptual framework of this linear theory is inadequate for the formulation 
of an exact theory. The appropriate conceptual framework was introduced 
in [5] and an exact nonlinear theory was proposed. The basic concept intro- 
duced in [5] is that of the material manifold which captures those properties 
of a crystalline solid which are intrinsic to it, being independent of the way 
it is embedded in space. The dynamics of the elastic body is then described 
by a mapping from space-time into this material manifold. 

While an exact theory of crystalline solids with a continuous distribution 
of dislocations was formulated in [5], the theory was left undeveloped up 
to the present time. The aim of the present work is to develop the theory 
and derive results which may be brought into contact with experimental 
data. Our focus in the present paper is on the static case, where we have 
a crystalline solid with a uniform distribution of elementary dislocations in 
equilibrium in free space. As a basic example we study, in the continuum 
limit, a two-dimensional crystalline solid with a uniform distribution of edge 
dislocations, one of the two elementary types of dislocations. In this case the 
material manifold is the affine group of the real line, the hyperbolic plane. 
We then study, in the continuum limit, a three-dimensional crystalline solid 
with a uniform distribution of screw dislocations, the other elementary type 
of dislocation. In this case the material manifold is the Heisenberg group. 

The purpose of the present work is to introduce to the mathematics and 
theoretical physics community a field where beautiful differential geometric 
structures, in particular Lie groups, form the basis of a classical physical 
theory. Moreover, the laws of the theory form a nonlinear system of varia- 
tional partial differential equations, which in the static case is elliptic and 
in the dynamical case is of hyperbolic type. Both cases constitute a worthy 
challenge for the geometric analyst. 
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This paper is organized as follows. In Part I, we give a general introduc- 
tion to the theory of crystalline solids containing an arbitrary distribution of 
dislocations, based on the work of Christodoulou [4], [5]. We introduce the 
basic concepts of material manifold and crystalline structure. In the present 
work, the fundamental notion is the canonical form, which is used to define 
the dislocation density. We illustrate the theory by giving two basic exam- 
ples of solids with a uniform distribution of dislocations. In order to state 
the laws governing the dynamics, the thermodynamic space is introduced 
and the energy function defined. The dynamics is described by a mapping 
from space-time into the material manifold. The equivalence relations for 
crystalline structures and for the mechanical properties of a solid are dis- 
cussed to give the proper physical interpretation of the theory. Finally, the 
Eulerian picture is given including the non-relativistic limit. In the Eulerian 
picture the material manifold is eliminated. 

In Part II, we focus our attention on the static case. We derive the 
boundary value problem from an action principle and give the Legendre- 
Hadamard conditions for the energy function. We then consider the two 
examples of uniform distributions of elementary dislocations and motivate 
the choice of our model energy function. We separately discuss the special 
cases of uniform distributions of edge and screw dislocations in two and 
three dimensions, respectively. In concluding the second part, we derive the 
scaling properties of the theory. 

Part III is devoted to the analysis of equilibrium configurations of a crys- 
talline solid with a uniform distribution of elementary dislocations in two 
and three dimensions. We solve the problem in two dimensions and give the 
method for the solution of the anisotropic problem in the case of a uniform 
distribution of screw dislocations in three dimensions. 
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Part I The General Setting 

1 The Material Manifold 

Let M be an oriented n-dimensional differentiable manifold, called the ma- 
terial manifold. It describes a material together with those of its properties 
which are intrinsic to it being independent of the way in which it is to be 
embedded in physical space. A point y G M represents a particle. We denote 
by X{M) the C°°-vectorfields on the material manifold Af. Let 

^' X ^ eyiX)=X{y) 

be the evaluation map at a point y G J\f. 

Definition I.l. A crystalline structure on N is a distinguished linear 
subspace V of such that the evaluation map ey restricted to V is an 

isomorphism for each y G M. 

Remark I.l. The orientation of M induces an orientation in V such that 
€y is orientation preserving. N admits a crystalline structure if and only if 
N is parallelizahle, see [5] . 

We introduce a 1-form v od. N with values in V defined by 

(2) Uy{Yy) = ey\Yy) ^Yy TyM . 

In the case that (A/", V) is a Lie group with corresponding Lie algebra, v is 
the Maurer-Cartan form. 

Remark 1.2. The most fundamental object is the 1-form v. We may in 
fact replace V by an abstract real vectorspace W of the same dimension, n, 
as the manifold J\f. Then we define a canonical form i' to be a W -valued 
1 -form on N such that 

Vy := v\^^^ -.TyM^W 

is an isomorphism for each y G A/". Given an element v G W , we may then 
define a vectorfield Y^ on M by 

Y,{y) = Uy\v) -.yyeM 

We then define the crystalline structure V by 

V = {X^ -.v (^W} . 

Then the canonical form corresponds to the 1-form v defined in 
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We say that the crystalhne structure V on is complete if each X £ V 
is a complete vectorfield on M. Then each element of V generates a 1- 
parameter group of diffeomorphisms of M, which represents (in the con- 
tinuum limit) a group of translations of the crystal lattice with parameter 
proportional to the number of atoms traversed. 

A complete crystalline structure V on the material manifold Af defines an 
exponential map 

exp : TV X V TV 

as follows. Let exp(y, X) be the point in TV that is at parameter value 1 
from y along the integral curve of X initiating at y. For each y G M, let 

expy : V TV 

X ^ expy(X) = exp(y,X) . 

We have 

exp^(O) = y , de^py{0) = ey . 

Thus dexpy{0) is an isomorphism for each y € TV. By the implicit function 
theorem it follows that, for each y G TV, there is a neighborhood Uy of the 
zero vector in V such that exp^^ restricted to Uy is a diffeomorphism onto its 
image in TV. 

Now choose a totally antisymmetric n-linear form a; on V which is positive 
when evaluated on a positive basis. The n-form to defines a volume form 
dfiuj, called mass form on TV by: 

(o\ dfM^{Yi^y,...,Yn,y) = U}{€-^{Yi^y),...,ey\Yn,y)) 

■.yyeAf,yYi^y,...,Yn,yeTyAf. 
The volume assigned by dfiuj to a domain V C TV, 




is the rest mass of P. 

Definition 1.2. Given a crystalline structure V on M we can define a map- 
ping 

A : AA^£(VA V,V) 

by: 

(4) A{y){X,Y) = ey'{[X,Y]{y))eV , yyeAf,X,YeV. 
We call A dislocation density. 



Suppose X,Y £ V are complete and generate 1-parameter groups of dif- 
feomorphisms ^'i, t G M from TV to TV. For y £ M 



D. CHRISTODOULOU AND I. KAELIN 



5 



coincides with (t, s) i— t- H^^, where is the 1-parameter group of diffeomor- 
phisms of M generated by A{y)(X,Y). 

The dislocation density is a concept that arises in the continuum limit of 
a distribution of elementary dislocations in a crystal lattice. An elementary 
dislocation has the property that, if we start at an atom in the crystal lattice 
and move according to one group of lattice transformations k atoms in one 
direction, then according to a different group / atoms in a second direction, 
according to the first —k atoms and finally — / atoms in the second directions, 
then we arrive at a different atom than we started from, but which, provided 
the circuit encloses a single elementary dislocation, is reached at in a single 
step corresponding to a lattice translation. This step is called Burger's 
vector. 

If A is constant on Af then, for all X,Y £V, there is a Z G V such that 

[X,Y] = Z, 

corresponding to a uniform distribution of elementary dislocations of the 
same kind, see below. Thus V constitutes in this case a Lie algebra, i.e. a 
vectorspace V over M with a bracket 

[., .] : V A V ^ V 

satisfying the Jacobi identity. 

By the fundamental theorems of Lie group theory, upon choosing an iden- 
tity element e £ Af, the material manifold Af can then be given the structure 
of a Lie group such that V is the space of vectorfields on Af which generate 
the right action of the group on itself. V is then the space of vectorfields 
on M which are left invariant, i.e. invariant under left group multiplication. 
The dual space V* is then the space of left invariant 1-forms on J\f. 

Let us consider diy, a 2-form on J\f : for any pair of vectorfields X, Y on 
M we have 

(5) du{X, Y) = X{u{Y)) - Y{u{X)) - u{[X, Y]) 

by the formula for the exterior derivative of a 1-form. In particular, this 
holds for X,Y G V. Now for X G V we have 

u{X){y) = Uy{Xy) = ey\Xy) = X , yyeM,XeV, 

a constant V-valued function on M. Similarly with X replaced by Y. There- 
fore, from ([5]) we have 

(6) diy{X,Y) = -u{[X,Y]) ,yX,YeV. 
On the other hand, 

(7) Aiy){X,Y) = e-\[X,Y]iy)) = z.,([X, y](y)) , Vy G JV,X,Y G V. 
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Thus, comparing ^ and ([7D, we obtain 
(8) {duiX,Y)) (y) = -Uyi[X,Y]{y)) = -I 

Let us then define the V-valued 2-form X on M by: 



My){x,Y). 



>^y{Yi,y,Y2,y) = A{y) {ey\Yi,y),eyHY2,y)) : yY,,y,Y2,y G T.AA 



at any point y G Af. We conclude from ([8]) that 



di' = —A . 



Let r be a closed curve in Af and let S be any surface spanning F, i.e. 9S = T. 
We finally conclude 



The right-hand side of ([9]) is the sum of all Burger's vectors enclosed by the 
curve r (or threading S). 

1.1 Uniform Dislocation Distributions and Lie Groups 

At the atomic level, two kinds of elementary dislocations are found. They 
are called edge and screw dislocations. In the following, we consider these 
two types of dislocations taking them as our model cases. For a uniform 
distribution of these two types of elementary dislocations, we determine 
the corresponding Lie groups, the affine group and the Heisenberg group, 
respectively. For a detailed description see [8] . 

1.1.1 Edge Dislocations and the AfRne Group 

The most basic type of a dislocation in a 2-dimensional crystal lattice is 
an edge dislocation. It appears in a 2-dimensional lattice in which an extra 
half-line of atoms has been inserted along the positive 1st axis. A circuit 
of translations in the directions of the 1st and 2nd axis, alternately, which 
encloses the origin, ends at an atom which is reached in a single step by a 
translation in the direction of the 2nd axis. On the other hand, circuits not 
enclosing the origin close. Mathematically, this phenomenon is represented 
by the commutation relation [Ei, E2] = E2, where E\,E2 are the vectorfields 
along the coordinate axis. 

We want to show that a uniform distribution of edge dislocations in a 
2-dimensional lattice gives rise in the continuum limit to the afRne group. 
This group is characterized by transformations of the real line of the form 



(9) 




,,,1 9 
X I— )■ e'' X + y 
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where {y^,y^) G are two parameters. The subgroups of the affine group 
are t ^ e^^t (multiphcation) and s ^ s + y"^ (translation) . 

We have as the group manifold equipped with the following multipli- 
cation 

{y\y^){y\f) = {y'+y\y'' + ey'f). 
X = A and Y = ey"^ 
generate the right action with 

' dy^ dy'^ dy"^ dy^ dy^ 

The Lie algebra of the affine group is thus generated by the vectorfields 
X, Y, which satisfy the commutation relation 

[X,Y]=Y. 

Therefore, the affine group is the material manifold J\f endowed with the 
crystalline structure. 

If we take {Ei = X,E2 = Y} as a basis of V, we have the following dual 
basis {oj^jOj^} for V* 

uj^=dy^ , uj^ = e-y'dy^. 
The corresponding left invariant metric 

(10) n= (u;l)2 + (u;2)2 ^ ^^^1)2 ^ ^^2y^ (^^2)2 

on M makes M the hyperbolic plane. 

The dislocation density A is 

X{Ei,E2){y) = Aiy)iEi,E2) = e^' {[Ei,E2]iy)) = {eyiE^)) = E^ . 
Since wi A a;2 = e"^^ [dy"^ A dy^) , it follows 

A = e-y^ [dy^ A dy2) E2 , 

and therefore 

f \= f e^y^ {dy' A V) E2 = A{J:)E2 , 

where A{T,) is the area of the surface S, a domain in Af. This makes sense 
since the sum of the Burger vectors associated to a domain in a uniform 
distribution of edge dislocations should be proportional to the area of the 
domain. 
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1.1.2 Screw Dislocations and the Heisenberg Group 

The second kind of elementary dislocation is called a screw dislocation. It 
appears in a 3-dimensional lattice in the following way. A circuit of transla- 
tions along the direction of the 1st and 2nd axis, alternately, which encloses 
the 3rd, ends at an atom which is reached at a single step by a translation 
in the direction of the 3rd axis, while circuits not enclosing the 3rd axis 
close. Mathematically, this phenomenon is represented by the commutation 
relations [Ei,E2] = E3, [Ei,E3] = [£"2,^3] = 0, where Ei,E2,E3 are the 
vectorfields along the coordinate axis. 

We want to show that a uniform distribution of dislocations of the screw 
type give rise in the continuum limit to the Heisenberg group. This group 
is represented as a group of unitary transformations on the space of square 
integrable complex valued functions ^ on M as follows 

L2(M,C) ^ L2(E,C) 

where {y^,y'^,y^) £ are three parameters. The subgroups of the Heisen- 
berg group are 1 1— t- '^{x -\- 1) (translation in position), s 1— )• e*'*^^(x) (trans- 
lation in momentum), and u 1— )• e*"^'(3;) (multiplication by a phase). 

We have as the group manifold equipped with the following multipli- 
cation 

iy\y^y^){y\f,f) = {y^ + y\y^ + f,y^ + f + y^f) ■ 

9 <9 id . ^ d 

generate the right action with 



d d id 
dy^ ' dy"^ ^ dy^ 



^ = Z , [X,Z] = [Y,Z] = 0. 



dy 



The Heisenberg group is thus generated by the vectorfields Ei = X, £'2 = 
Y, £3 = Z, which fulfill the commutation relations 

[Ei,E2] = £3 , [£i, £3] = , [E2, £3] = , 

and generate the right multiplication. The linear span of {Ei, £21 £3) forms 
a Lie algebra (crystalline structure) corresponding to a uniform distribution 
of screw dislocations in a three-dimensional crystal lattice. Therefore, we 
can associate the Heisenberg group, which is the group corresponding to the 
crystalline structure, with the material manifold M. 
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If we take {Ei, £^2,-5/3} as a basis of V, we have the following dual basis 
{iu\ 0^^,00'-^} for V* 

= dy^ , a;2 = dy^ , uj^ = dy^ - y^dy^ , 

and the corresponding metric 

(11) n= {u'f + {U^f + (^3)2 ^ (^^1)2 + (^^2)2 + (^y3 _ ^1^^2)2 ^ 

which is a Bianchi type VII metric. The manifold M endowed with this 
metric is a homogeneous space. 

Here, the dislocation density A turns out to be 

\[Ei,E2)[y) = E^,X{EuE^){y) = \{E2,E^){y) = Q, 
and therefore 

\ = {uJl^ UJ2) E-s = {dy^ A dy"^) E3 . 
The integral of A over a surface S in A/" is 

/ A = / {dy^ A dy^) Es = A(nS)i?3 , 

where 11 is the projection map of the line bundle of the homogeneous space 
pip over M? with the standard metric on the base (the curvature of the 
bundle being —dy^ A dy'^). 

2 The Thermodynamic State Space 

Consider the space 5'2'(V) of inner products on the crystalline structure V. 
The thermodynamic state space is defined as the product 

and its elements are (7,0"), where 7 G 5^(V) is the configuration and 
a £ M"*" is the entropy per unit mass. 

Each 7 G 5'^(V) defines a totally antisymmetric n-linear form on the 
crystalline structure V by the condition that if (Ei, . . . , En) is a positive 
basis for V, orthonormal with respect to 7, i.e. -yAB ■= i{Ea,Eb) = 5ab-, 
then 

u;y{Ei, . . . , En) = 1. 
It follows that there is a positive function V on S'^(V) such that 

uj^ = ^^(7)1^ . 

The positive real number V{j) is the volume per unit mass corresponding 
to the configuration 7. 
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The thermodynamic state function k is a real-valued function on the 
thermodynamic state space ^^(V) x M"^. The Lagrangian which determines 
the dynamics will be defined through this function. 

The thermodynamic stress corresponding to a thermodynamic state 
(7,0") is the element 7r(7,cj) of (5'2(V))* defined by 

— = --^(7, a)V{^) . 

The temperature corresponding to a thermodynamic state (7, a) is the 
real number 'd{j, a) given by 

with the requirement that 'd{'y, a) is positive and tends to zero as a tends 
to zero. 



3 The Dynamics 

In the general theory of relativity the space-time manifold is an oriented 
{n + l)-dimensional differentiable manifold Ai endowed with a Lorentzian 
metric g, that is a continuous assignment of a symmetric bilinear form 
of index 1 in T^M ,\/x £ M. The Lorentzian metric divides TxA4 into three 
different subsets Ix, Nx, Sx, the set of timelike, null and spacelike vectors 
at X respectively, according to whether gx restricted to the corresponding 
subset is negative, zero or positive. The subset Nx is a double cone called 
the null cone at x. The subset Ix is the interior of this cone, an open set of 
two components, the future and past component. Sx is the subset outside 
the null cone, a connected open set for n > 1. A curve 7 is called causal 
if its tangent vector belongs to Ix U Nx, Vx G 7, and it is called timelike if 
its tangent vector belongs to Ix, G 7. We assume that {M.,g) is time 
oriented, that is a continuous choice of future component Ix of Ix can and 
has been made Vx G Al. This choice determines the future component Nx 
of Nx at each x £ M. A timelike or causal curve is then future or past 
directed, according to which component its tangent vector belongs to. A 
hypersurface Ti is called spacelike if at each x £ Ti the restriction of gx to 
TxTi is positive definite. A spacelike hypersurface in Ai is called a Cauchy 
hypersurface if each causal curve in A4 intersects "H exactly once. We 
assume that {M.,g) possesses such a Cauchy hypersurface [3j. 

The motion of the material continuum is described by a mapping / from 
the space-time manifold A4 into the material manifold M, 

(12) f:M^M. 
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This mapping specifies which material particle is at a given event in space- 
time. It is subject to the following conditions: 

i) The mapping /, restricted to a Cauchy hypersurface 7i, /|^, is one to 
one. 

ii) The differential of the mapping /, df{x), has a 1-dimensional kernel 
contained in Ij;, Vx S A^. 

Then, for each y £ fi-M) C A/", f~^{y) is a timelike curve in Ai. The 
material velocity u is the future directed unit tangent vectorfield of the 
timelike curve f~^: 

span(u^) = ker(4f(x)) = T^f'^iy) , /(x) = y, g{ux,Ux) = -l,yx £ M . 

The simultaneous space at x is the orthogonal complement of the linear 
span of Ux- 



Note that gj]^ = gx\-£ is positive definite. The restriction of the differential 
df{x) to T,x is an isomorphism of S^. onto TyM, where y = f{x). 

iii) The isomorphism df{x)\j,^ is orientation preserving. 

The equations of motion, a second order system of partial differential 
equations for the mapping /, are derived from a Lagrangian L, a function 
on M which is constructed from /. The action ^ in a domain T) d M is 
the integral 



where d^g is the volume form of {M , g) . 

Any mapping / fulfilling the three requirements stated above, defines an 
orientation preserving isomorphism jf^x of V onto T,x by 



T,x = (span(«2;)) 



± 




(13) 



jf,x = {df{x)\j.J 



-1 



Define 



(14) j}^x9x = 7 , 

7 G ^^(V) and depends only on / and x. Note that 



1 




is the rest mass density at x. 
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To take into account thermal effects in the Lagrangian picture we invoke 
the adiabatic condition which states that the entropy per unit mass of any 
element of the material remains unchanged. This allows us to consider the 
entropy per unit mass cr as a given function on the material manifold M. 

The Lagrangian function L on 7W is defined by: 

(15) L{x) = K {jl,,g,,a{f{x))) ,yxeM, 

where n is the thermodynamic state function on S'^(V) x M"*". Note that 
L{x) depends on g only through and not on derivatives of g. 

The energy-momentum-stress tensor Tx at x is an element of the 
dual space {S2{TxJ^))* defined by 

(16) g(L(x)d;.,(x)) ^_l 

Ogx ^ 

Prom [5] we have the following 

Proposition I.l. We can write the energy-momentum-stress tensor at x as 
follows: 

(17) Tx = p{x)ux 0Ux + Sx 

with the mass- energy density p (since we are in the relativistic framework 
this includes the rest mass energy) given by 

(18) pix) = Kijl^gx,a{f{x))). 
The stress tensor is given by 

Sx{9x) = ■^{ff,x9x,crifix))){j}^x9x) ■ 

Since jfx9x = j} x9^x> can be viewed as an element of (52 (S^))*- 

Remark 1.3. Hence we see from Proposition \I.1\ that in the case of pure 
continuum mechanics (in the absence of electromagnetic fields) the energy 
per unit mass is e = nV. 

Define the principal pressures pi, i = l,...,n as the eigenvalues of 
S\,[, relative to ^l^, where the subscript bb means lowering of the indices 
with respect to g. Note that the principal pressures pi : i = 1, . . . ,n can 
equivalently be described as the eigenvalues of tt;,;, relative to 7. 

The positivity condition on Tx requires that 

H+ ^ TxM 

(19) ^ -Tlfv" 
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maps into Ix ■ This is equivalent to the condition that \pi\ < p : i = 
l,...,n. Here, in the case of crystahine soUds, we assume the stronger 
condition that the range of (|19p hes in which is in turn equivalent to 
\pi\ < p : i = 1, . . . ,n. 



3.1 Variation of the Lagrangian 

Let V = df{x) G U(^x,y) ^ C{TxM.,TyM), y = f{x), where U(^x,y) is the open 
subset of the linear space C{TxAi.,TyJ\f) consisting of those v G C{TxAi,TyJ\f) 
which verify the two conditions 

1) keru is a time like line in T^-TW, 

2) with Tlx the (72,-orthogonal complement of ker-y in TxM, the isomor- 
phism 

is orientation-preserving. 

These above conditions 1) and 2) correspond to the conditions ii) and iii), 
respectively. Note that i) is not a local condition so it does not reduce to a 
condition on v. 

Since the isomorphism 

as defined in ([T3|), depends only on we may write 

i{v) = jf,x , 

so the configuration jfx corresponding to the isomorphism v £ ^{x.y) is 
given by the positive quadratic form on V 

(20) ^{v)=e{v)g\^^ . 
Now V I— 7- 7(u) is a mapping of 

(21) ^■■= U ^(-.?/)' 

{x,y)<^My.M 

into 5^(V). The bundle 5, a bundle over M. x A/", is the bundle over which 
the Lagrangian is defined. 

The mapping v 1— )• ^{v) is described by the functions ^ab{v) on where 

(22) -iAB{v)=i{v)[EA,EB) , 

where Ea : A = 1, . . . , n is a basis of V such that . . . , En) = 1. 
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By the properties of v we have a positive basis (Xj^^x) : A = 1, . . . , n) of 
T,x defined by: 

V ■ Xa{x) = EA{y) : A = l,...,n, 
where EA{y) = ey{EA) G TyN{A = 1, . . . , n), and thus from ([20]) and ([22]) 

lAB = 5Ie, {Xa.Xb) ■ 

Note that (n, Xi, . . . , is a frame field for M. Let v G C{TxM,TyM) 
be a variation of v S U^^^^yy To describe ■)) we must give v ■ u £ TyM and 

v-'u(x) = v^EAiy), 
v-Xa{x) = v^EBiy). 

This is because {EA{y) : A = 1, . . . , n) is a basis for TyM . So the {vq : A = 
1, . . . ,n;v^ : A,B = 1, . . . , n) can be thought of as the the components of 

V. 

In view of and the Lagrangian function L is L{v) = p {'y{v),a{y)), 
where p is the relativistic energy-density which includes the rest mass con- 
tribution, a function of a thermodynamic state {'y{v),a{y)) € 5'^(V) x M+. 
L is a function on the bundle B (|2ip . The Lagrangian form, a top degree 
form on A^, is Ldpg. 

The Lagrangian L depends on v through the configuration 7, L^jab), 
and the first variation reads 

■ dp. 

L = ^ 7AB , 

where the first variation of ^ab is given by 

lAB = -VaIBC - VbIAC 

(see [4]). To formulate the hyperbolicity condition (see below), we need to 
consider the second variation of L with respect to v, 

d'^L 

(23) L = ■ {v,v) = h{v,v) , 
where, in general, 

(24) h{v, v) = h%v^v^ + 2h%v^v^ + h'^EvcvE ■ 

For L = L{'^ab)-, using the formula for the second variation of ")ab^ 

lAB = 2 {lAClBDv'^V^ + ICDVaVb + IACVb^d + IBCVaVd) 
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(see [1]), we obtain 

" dp .. d'^p 

= ^ lAB + lABlCD 

OlAB OJabOJCD 

- {-fAClBDVQVQ + IACVbVd + IBCVa^D + ICDVaVb) 



AB 



q2 

+ 75 7^ {vflEB + VbIEa) {vclFD + V^IFC) 

ojabojcd 

dp An dp An 

2^ IacIbdVq Vq + 2- ^abVcVd 

O^CD OJCD 

, „ xc\ -A-B 



+2 -TT^lBE^'i + TT^IAeS'b VcVd 
\OJCE O-fDE J 

d^P -A-B 

+4^ t: lAElBFVcVr) ■ 

d^CEd^DF 

Comparing coefficients witli (j24p using (j23p . we see tliat: 

h'L = 2/^jAC7BD, h<i%=0, 
OlCD 



11 u A d'^P 

I^AB = j: lAElBF 

d-fCEd-iDF 



.CD 



+2 (^^^AB + J^IBeS^ + TT^TAE^g ) , 
\djCD d-fCE djDE J 

wliere tlie last expression will be made use of in the derivation of the 
Legendre-Hadamard conditions in Part [III 



3.2 Hyperbolicity and Characteristic Speeds 



Set 

n — hOO n/fCD _ uCD 

'-'AB — -n-AB ) ^^AB — I^AB ■ 

Then 

h{v,v) = -GabVqVo + M'^EvcVd ■ 
The hyperbolicity condition in general is that there is a pair (T, 9) E T^M x 
T*M with 9 ■ T > such that h is negative- definite on 

(25) {i -.v = 9^Y,Y £ TyN} 
and positive- definite on 

(26) {i):{; = K®y,K-r = 0,yG TyM} 

(see [4]). We stipulate here that the above conditions hold with (T, 9) defined 
by the rest frame of the material at x, i.e. T = Ux, = 
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For V in (|25|) we have: 

V • u = Y , V • Xa = , 

and the condition that h is negative definite on (I25p is the condition that 

(27) GabI^^^^ > : Vy / , 
i.e. that G is positive-definite. 

Since 

Gab = T^AB + PlAB 
(we are raising and lowering indices with respect to 7ab) and the principal 
pressures pi, . . . ,pn are the eigenvalues of ttab with respect to ^ab , this 
condition is: 

minpi > -p, 

i 

which of course follows from maxj \pi\ < p. For v in (|26p we have: 

V • u = , V ■ Xa = k-aY , 

where ka = k ■ Xa, i-e. 

v^ = , v^ = kaY^, 
and the condition that h is positive definite on (j26p is the condition that 

(28) M^I^kcKdY^Y^ > : Vk / 0, Vy / . 
Set 

Then condition (j28|) is that HAB{f^) is positive-definite for all k 7^ 0. 

If the first condition (j27p is satisfied, the second condition ()28p becomes 
the condition that for k 7^ the eigenvalues Ai(k), . . . , A„(k) of Hab{i^) 
with respect to Gab are all positive. Note that Hab{i^) is homogeneous of 
degree 2 in k, hence so are the Ai(k), . . . , A,i(k). 

Now the characteristic matrix is 

Xab{€) = ^AB^^^iy ■ 

We use the basis {u, Xi, . . . , X„) for T^J^- We denote by w the frequency 
and by ka, A = 1, . . . ,n, the wave number components, i.e. 

■=Co-=^-u , ka:=U-=^-Xa- 

Then 
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so 

n 

det X = det G • JJ (Ai(K) - w^) . 

i=l 

We see that the second condition is equivalent to the 2n roots ^^J\i{^^), 
i : l,...,n, of the characteristic polynomial detx as a polynomial in u 
being real. The characteristic speeds are 

(29) r/i = , where |k| = J (t"!)^'^ kai^b ■ 

\k\ * 

Note that the ryj are homogeneous of degree in k. The causality condition 
is that the inner characteristic core in T*M contains the null cone of g in 
T*M. This reads (in units c = 1) 

m<i , Vz = 1, . . . ,n. 

4 Properties of the Energy per unit Mass 
Postulate I.l. We stipulate that the energy per unit mass e = pV has a 

o 

strict minimum at a certain inner product 7- 

Remark 1.4. The above postulate has an intuitive physical interpretation. 
Note that the set of inner products 7 is an open positive cone in the linear 
space of quadratic forms. For large 7 (large expansion) and also for 7 near 
the boundary (large compression) the energy e is physically expected to blow 
up, so we can restrict ourselves to a compact set of inner products, where e 
necessarily attains a minimum. 

o o 

We choose {Ei, ...,£'„) to be an orthonormal basis relative to 7, so 7yi_B= 
6ab- This is compatible with the previous condition uj{Ei, . . . , E^) = 1 for 
a suitable volume form u on V. This choice of uj corresponds to a choice of 

o 

unit of mass so that the mass density associated to the configuration 7 is 

o o . 

equal to 1. 7 defines a metric n on A/ by 

n n 

A,B=l A=l 
where {ui^ , . . . , w") is the dual basis to {Ei, . . . , En). If V is a Lie algebra, 

o 

so that A/ is a Lie group, then n is a left-invariant metric, i.e. invariant 
under the actions of M on itself by left multiplications by elements of N . 
Thus (AA, n) is a homogeneous Riemannian manifold (which in general is not 
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isotropic) . The Riemannian manifold {Af, n) has curvature except when V 
is Abehan, so there are no dislocations. 



Note that djio^ = Vdet nd^y = dfi^j, therefore we have for any domain 
in Af 



M{n) = 




n n 



that is, in our choice of units the mass of Q is equal to its volume with 

o 

respect to n. 

4.1 The Isotropic Case 

o 

Let us define the orthogonal group corresponding to 7, 

Oo = jo G C{V,V) -.7 {OX,OY) =7 {X,Y),yX,Y G v} . 
Oo acts on S'2~(V) in the following way: 7 1— t- O7, where 

(O7) {X,Y)=j{OX,OY). 

000 o 

Since 7 1— )• O 7 =7, the symmetric bilinear form 7 is a fixed point of the 
action of Oo on S'^(V). 

If the energy density 6(7) is invariant under Oo we are in the case of 
isotropic elasticity. Then 6(7) depends only on the eigenvalues Ai, . . . , A„ 

o 

of 7 relative to 7, 

e(7) = e(Ai,...,A„), 

where e(Ai, . . . , A^) is totally symmetric in its arguments. This is the sim- 
plest case of an energy density. 

Consider 

q = f n . 

This is a 2-covariant symmetric tensorfield on M, i.e. at each x & Ai is a 
quadratic form in T^M. The vector Ux belongs to the null space of Qx and 
the restriction QxIy:^ positive-definite. 

Let Ai, . . . , An be the eigenvalues of qx\-£^ relative to We then have 

Proposition 1.2. The eigenvalues Ai, . . . , A.„ are the inverses of the eigen- 

o 

values Ai, . . . , A„ of j = jj^gx relative to 7- In particular, 

VAi-----An= ^ ^ =- = N. 

• • • • • A„ V 
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Therefore, in the variational principle, the crystalline structure V on M is 

o 

eliminated in favor of the Riemannian metric n. 
5 Equivalences 

5.1 Equivalence of Crystalline Structures 

Definition 1.3. Two crystalline structures V and V' on J\f are said to be 
equivalent if there is a diffeomorphism ip of M onto itself such that -0* , the 
push-forward ofip, induces an isomorphism ofV onto V'. 

Let z/ be the canonical 1-form ^ associated to V and z^' the one associated 
to V'. We have 



We may take the above as the definition of pullback for V- valued 1-forms: 
It then follows that: 



with the natural extension of the above definition of pullback to V-valued 
2-forms, i.e. 



(30) 



U'idi: ■ Yy) 



y G V : G TyN , 
V'^yeV : d^'-Yy^T^^y)M . 



(31) 



A = V*A' 



\' {dip ■ Xy.dilj - Yy) = ilj^{\{Xy,Yy)) I VX, F G V , Vy G . 



In fact, we have for all X, 1" G V and y G N: 



A'(V'*X,V'*r)(V^(y)) = i/'([V^,X,V^,y](V'(2/))) 



= z/'(V*[X,y](V^(y))) 

= v'{d^-[X,Y]{y)) 

= ^.{u{[X,Y]{y))) 

= V*(A(X,y)(y)) , 



where we have made use of the fact that [V'*X, 4'*Y] = i^^^lX, Y]. 



5.2 



Equivalence of Mechanical Properties 



We investigate the question of the equivalence of the mechanical properties 
of a solid. A certain solid phase of a certain substance is described by an 
equation of state. 
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Let ^ : V — 7- V be a linear isomorphism. The group of ah such A 
is homomorphic to GL„(]R). Then A* : S'^(V') — )■ 5'^(V) is the induced 
isomorphism defined by 7 = A* 7', where 

-/{X,Y) =-f'{AX,AY) ,yX,YeV. 

The corresponding energy functions on /^^(V) and S'^(V') are denoted by e 
and e' , respectively. 

Definition 1.4. Two energy functions eonV and e' on V are equivalent if 
there exists a linear isomorphism A : V ^ V such that e'(7') = 6(7), where 
-f = A*i, for all 7' G S^{V). 

We illustrate the above definition of mechanical equivalence by giving 
an example in the isotropic case. Let 7 G 5'^(V) be given, and define 
M G £(V,V) by 

7(MX,y) = 7(x,y) ,vx,yGV. 

Similarly, given an isomorphism A as above (i.e. 7 = A*^'), we define M' G 
£(V', V) from the corresponding 7' G 5'^(V') as 

o 

7' (M'x', y') = 7'(x', y') , vx', y G V . 

o ° 

In the above, 7= A* 7'. Then M' = AMA'^. For, given any X\Y' G V 
let X, y G V be X = A-^X', Y = A'^Y'. Then 

7'(M'x',y') = 7'(x',y') = 7(^,>') =7 (Mx,y) 

o o 

= 7' (AMX, yiy) =7' {AMA~^X', Y') . 

Therefore, the eigenvalues A'^, . . . , of M' coincide with the eigenvalues 
Ai, . . . , A„ of M. In the isotropic case, e'{'y') = 6(7) means 

e'(Ai, . . . , A^) = e(Ai,...,A„,), 

where both sides of the above equation are symmetric functions, so the 
energy functions coincide, i.e. e' = e. 

The equivalence of the energy functions, however, does not fully capture 
the equivalence of having two solids of the same substance in the same phase. 
What is required in addition is to have the same equilibrium mass density 
of infinitesimal portions. In fact, it is the triplet (V,a;,e) which defines a 
solid with its mechanical properties. Two solids of the same substance in 
the same phase, but with possibly different dislocation structures, to the 
extent that they can be described by the same manifold M (which is true if 
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they are diffeomorphic) are defined by the triplets {V,u),e) and (V',a;',e'), 
respectively. Additionally, there is an isomorphism ^4 : V — >■ V, such that 

io = A* J , 

i.e. a;(Xi, . . . , X„) = ijj\AXx, . . . , AX^) : VXi, . . . , X„ G V as well as 

(32) e'(7') = e(7) , 7 = ^ V : V7' G 52+(V') 

in accordance with Definition II. 4[ Xhus if {E\ , . . . , En ) is a positive basis 

o 

for V which is orthonormal relative to 7, then 

LoJyE\, . . . , En) = 1 , 

while uj{Ei, . . . , En) = fJ-o- On the other hand, with E'- = AEi : i = 1, . . . ,n, 

o 

{E[, . . . , E'^) is orthonormal relative to 7', thus 

Uo {E[,... ,E'n) = 1 , 
7' 

while uj'{E[, . . . , E'^) = uj{Ei, . . . , En) = Ho- Consequently, the mass density 

o 

corresponding to 7' relative to the triplet {V',uj', e') is the same as the mass 

o 

density corresponding to 7 relative to the triplet (V,a;,e). 

Remark 1.5. If we have two triplets iy,Lo,e) and (V',a;',e') which corre- 
spond to the same substance in the same phase and, additionally, the iso- 
morphism A : V ^ V is of the form A = where ^ is a diffeomorphism 
of N onto itself, then the dislocation structures are also equivalent. In this 
case, if f : M ^ N is a dynamical solution of the problem corresponding to 
(V,a;, e), then f':=ipof:J^^Misa dynamical solution of the problem 
corresponding to {V' ,e'). 

5.3 Similarity 

Let now V be fixed, so that 5'2'(V) is also fixed. Two different energy 
functions e and e' may be related as follows. There is a constant a > and 
an isomorphism V — ?• V, defined by 

X ^aX : VX G V . 

This defines an isomorphism ^^(V) — )• ^^(V) by 7 1— )• 7', where 

7'(X, Y) = 7(aX, aY) = a'^-f{X, Y) : VX, Y G V , 

i.e. 7 I— )• 7' = a^7. 

According to the above discussion of equivalence of mechanical properties, 
we must have 6(7) = e'{'y'), V7 G S'^(V), and, additionally u' = a^uj for the 
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same material in the same phase, for 

io'{Xi,...,Xn) = u;{aX, , . . . , X^) = a^^i, ■ ■ ■ , X^) . 

Then the two theories (primed and unprimed) represent the same material 
in the same phase. 

Moreover, since V is identical in the two theories, the dislocation struc- 
tures may seem identical. However, for a given domain Q CC N", which has 
the same crystalline structure V\q (the restriction to Q of the vectorfields 
in V), the primed theory actually assigns a physical dislocation density 
times that of the unprimed theory. This is due to the fact that ui' = a^uj 
and the observation that in any dimension, the elementary dislocations are 
co-dimension two objects. Thus their density refers to the two-dimensional 
measure of a cross-section. 

6 The Eulerian Picture 

As a prelude to formulating the Eulerian picture, we view the crystalline 
structure V as an abstract n-dimensional real vector space, as in Remark 
II. 2[ We then view the canonical form v as a 1-form on M with values in V 
such that i^lxyj^f is an isomorphism from TyM onto V for all y € M. Thus, 
given any v £ V, we have a tangent vector 

[i^It^^Y' ■V = Yy{v)eTyJ\f -^y^-^' 

that is a vectorfield Y{v) on J\f. Then the crystalline structure in the original 
sense is the space 

{Y{v) :veV} 

of vectorfields on J\f. 

To obtain the Eulerian description we must eliminate the material mani- 
fold M. A fundamental variable is the material velocity u, a future-directed 
timelike unit vectorfield on . This defines the distribution of local simul- 
taneous spaces 

(33) S = {S^ : X G A^} , 

where T,^ is the orthogonal complement of Ux in TxAi- 

We then need another entity defined on to play the role of the canonical 
form ly. Consider 

e = r^. 

This is a 1-form on M with values in V, viewed as an abstract n-dimensional 
vector space. 
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The 1-form ^ has the fohowing properties: 

1) ^^x = 0, 

2) for each x £ A4, is an isomorphism from Tix onto V, and 

3) Cu^ = 0. 

We thus introduce ab initio a 1-form on M with values in V possessing the 
above three properties, as another fundamental Eulerian variable besides u. 

By 2), given any v £ V, we have a tangent vector 
at each x £ M, that is we obtain a vectorfield X{v) whose value at each 



point belongs to the distribution (p3]l . i.e. it is orthogonal to u. 

A mapping A4 — )■ S'^(V), x i— )• 72;, is then defined by: 
(34) 

7x{vi,V2) = gxiXxivi),Xx{v2)) = g\j^^iXr,{vi),Xx{v2)) ■ yvi,V2£V. 
The last fundamental variable is the entropy s, a positive function on Ai. 

The volume form uj, the space of thermodynamic configurations «S'^(V), 
and the volume per unit mass V are defined as in Section [2l The ther- 
modynamic state space is 5*^(1^) x M"*" and the energy per unit mass e is 
a function on this space, as before. The thermodynamic stress vr at each 
(7,s) G S+{V) X M+, 7r(7,s) £ (52(V))*, is defined by 

as in Section [2j Also, the temperature 9 is defined, as before, by 

ae(7,s) 

Thus 

de = —-Vtt • d'y + Ods 
expresses the first law of thermodynamics in the present framework. 

The stress 5 is a 2-contravariant symmetric tensorfield on A^, an assign- 
ment of an element Sx £ {S2iX'x))* at each x £ M. This is defined as 
follows. Given any Qx £ S2{TxA4), we define ^x £ S2{V) as in (p4|) by: 

'yx{vi,V2) = gx{Xx{vi),Xx{v2)) = g\j^^ {Xx{vi),Xx{v2)) : yvi,V2 £ V. 

Then 

Sx{gx) = vr(7^,s(x))(7x) . 
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The mass-energy density p is then the positive function on Ai given by 

e(7^,s(x)) 

and the energy-momentum-stress tensor is defined according to ([T 



The Eulerian equations of motion are a first order system of partial dif- 
ferential equations consisting of 

(a) CuS, = 0, i.e. property 3) of ^, and 

(b) V-T = 0. 

Remark 1.6. In n space dimensions, dimA^ = n + 1, there are + n + 
1 dependent variables, the r? algebraically independent components of ^, 
the n algebraically independent components of u, and also s. The Eulerian 
equations are also n'^ + n + 1 in number, r? independent equations in (a) 
and n + 1 independent equations in (b) . 

For solutions of the Eulerian equations such that u, ^, and s are continu- 
ous, (b) implies the adiabatic condition on s: 

u{s) = 0. 

For the proof of this fact see [5] . 

Remark 1.7. The case of absence of dislocations is the case dS, = 0. 

Let {Ea : ^ = 1, . . . , n) be a basis for V. Given any 7 G 5'2'(V), we set 

lAB = ^{Ea,Eb) ■ 

According to the above, at each x ^ M there is a unique Xa S such that 

i-XA = EA ■.A = l,...,n. 
Thus Xa is a vector field on M belonging to S. Then 

S = tt^^Xa ® Xb , 
where tt^^ G S'^(V)* is the thermodynamic stress defined above. 

Fix an element v £ V and consider the vectorfield X C S such that 

C-X = v. 

Then 

= u{^-X) = {CuO-X + C-[u,X] , 
i.e. (according to 3)) ^ • [u, X] = 0. It follows that 
(35) U[u,X]=0, 
where H is the (7-orthogonal projection to S. 
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Proposition 1.3. Define := y. Then the mass current I = fiu satisfies 
the equation of continuity: 

(36) V • / = . 

Proof. V is endowed with a volume form lo and uj-y = Vlo. Choose a basis 
(£"1, . . . , En) for V such that uj{Ei, . . . , £'„) = 1. Then V = ^/detJ, and, by 



(37) JAB = i{Ea, Eb) = g{XA,XB) , 
hence 

(38) uiV) = ^V{^~yui^AB), 



and, using ([3 

uijAs) = u{g{XA,XB)) 

= g{VuXA, Xb) + g{XA,VuXB) 

= giVx^u, Xb) + g{XA, Vx^u) + g{[u, Xa],Xb) + g{XA, [u, Xb]) 
By (p5|) . the last two terms in the above equation vanish. Let us define k by 
(39) VxaU = k^Xb . 

[{Xa : a = 1, . . . , n) is a basis for Sj; at each point.] Then 
gC^XA^, Xb) = ka9{Ec, Eb) = ^cbka ■ 

Hence we obtain 

ui'jAB) = ICBKa + lAC^B ■ 
Substituting in ([38]) yields 



Then /_i = satisfies 

u(/i) + /itr K = , 
and since, from ()39p . tr k = V • u, we obtain 
(40) V • / = V(^m) = u{^jl) + ^iV ■u = 0, 

which establishes (1361). □ 
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6.1 The Non-relativistic Limit 

First of all, we restrict ourselves to the case where {A4,g) is the Minkowski 
space-time. Then we consider the non-relativistic limit, where Minkowski 
space-time is replaced by Galilean space-time. There, we have the hyper- 
planes T,t of absolute simultaneity, which are isometric to n-dimensional 
Euclidean space. Any family of parallel lines transversal to the repre- 
sents a Galilean frame, that is, a family of observers in uniform motion and 
at rest relative to each other. Any such family of parallel lines defines an 
isometry of the T,t onto each other. 

In terms of a Galilean frame and a rectangular coordinate system (x^, . . . , x") 
in Euclidean space, and with = ct (c: the speed of light in vacuum), the 
space-time velocity u = u^-^^ is represented in terms of the space velocity 



V = v'£, by 





S^^u^ = , u^ = g^^u^" , 
reads (note that u^, = g^uu'^i i-e. uq = —u^, while Ui = u^) 



(42) S'^-S'^— = 0, 

c 

(43) ^oo_^oi!^ ^ Q 



So, in the limit c — )• oo, we have 5*" = S^^ = and, therefore, 
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e = + e' , 

p = fj,e = p(? + e , e = fie' . 
We have for the mass current I'^ = fiu'^: 

/° = , = = u+ -^-^ + 0(c- 



Vl-|?;|2/c2 c2 2 c3 

In the non-relativistic hmit the equation of continuity V ^1^ = becomes 
the classical continuity equation: 

Similarly, for the energy- momentum-stress tensor T^a, we have 
TOO = p{uy + = pc^ + {e + pv') + 0(c-2) , 

c 



The equations of motion VxT'^ = reads 

dx^ dx^ dt dx^ 

dx^ dx^ dt dxi 

In order to obtain the energy equation, let us consider 

For the components of F'^ we have, using (|42|) . i.e. 5**^ = S^^v^c and S^^ 

^ = 1 12 / 2 + 2 , ,0 , o = + C"2 



(^|^|2 + ^ + 5^^- - i/ib|2- + 0(c"2) 

c c 2 c 
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Hence, 



is the energy equation 
d (\ 



QpO Qpi 



(45) ^U^H^ + ^)+^ 



d 



0. 



The non-relativistic Eulerian equations are 
(a) -± + C,^ = 0, 
rill. r) 

0, 



(bO) 

(bl) 

(b2) 



dt dx' 



dt 



d_ 

dt 



dx^ 

-Li\v\'^ + e I + 



1 



_d_ 

dx^ 



1 



0. 



Notice: 



• (bO) is the differential mass conservation law, a consequence of (a), 

• (bl) is the differential momentum conservation law, and 

• (b2) is the differential energy conservation law. 

Remark 1.8. (a) may be thought of as expressing a law of conservation of 
dislocations. In fact, let Cq be a closed curve on Sq. Let (pt be the flow in 
Galilean space-time generated by the vectorfield 



d 



dt dx^ 

Then (ptl^^ is a diffeomorphism from So onto Sj. Let Ct = (ptiCo). Then, 
according to (a), 



Ct Co 

The left-hand side corresponds, in the continuum limit, to minus the sum of 
the Burger's vectors of all the dislocation lines enclosed by Ct- 

Remark 1.9. For continuous solutions, i.e. v, s) continuous, (b2) is 
equivalent to the adiabatic condition 



modulo the other equations. When discontinuities such as shocks develop, 
this equivalence no longer holds. However, the Eulerian equations (a) -(b) 
still hold, but in a weak or integral sense. 



ds 
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6.2 Prom the Eulerian to the Lagrangian Picture 

In order to go from the Eulerian picture to the Lagrangian formulation we 
have to extract the canonical form u from ^. First note that 

(46) £„^ = 4^ 4>*tC = C, 

where (pt is the flow generated by u. 

Let Tihe a Cauchy hypersurface in M. We identify T-L with M and define 
f : M ^ J\f as follows: 

fix) = yen 

is the point at which the integral curve of u through intersects T-L. 



For Xx G TxM we have from (|i6]) 

Define v to be the V-valued 1-form induced by ^ on "H: 

We must show 
Proposition 1.4. 

ri^ = ^- 

Proof. Let Vp € TpA4, p G (pti'H) for some t G M. uniquely decomposes 
into 

Vp = Xu + Yp , where Yp G Tp(f)t{T~i) ■ 
Therefore, it suffices to show that {f*i^){up) = (,{up) = 0, which is obvious, 
and 

(/V) (Yp) = ^Yp) . 
Since is a diffeomorphism from T-L onto (ptiT-L), dcptiq) is an isomorphism 
from TqH onto Tp(j)t(T-l), 4>t{<l) = P- Therefore, there exists a unique Xq G 
TqH, such that 

d(l)t ■ Xq = Yp . 

Now, 

i{Yp) = mt ■ Xq) = {Xq) = i{Xq) = u{Xq) 

and 

(/*^) = (/*^) • = {(blf*!^) {Xq) = (/V) {Xq) = U{df . Xq) . 

But f\y^ = id, hence df • Xq = Xq because 

df\T,H = ^d. 

□ 
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Remark 1. 10. The Eulerian picture is the one more related to direct phys- 
ical experience and also the one which may serve as a basis for extending 
the theory beyond the domain of elasticity theory, where the dislocations are 
no longer anchored in the solid. We shall see in the next part, however, 
that the Lagrangian picture provides the suitable framework for the study of 
static problems. 
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1 Formulation of the Static Problem 



Let M = ft C M."" be the material manifold and Ai = -E" Euclidean space 
(n = 2, 3 being the cases of interest). In the static case, we adopt the ma- 
terial picture, that is, we consider one to one mappings (p from the material 
manifold into space, 

(47) "^'-^ ^ 

y ^ <p{y) = x. 

They correspond to mappings f ■ M x R — )• A/" such that 

f{x,t) = ^~\x) :Vt. 

In the material picture, the interchange of the roles of the domain and target 
space transforms a free boundary problem into a fixed boundary problem. 

We recall that the configuration 7 is an element of the space of inner 
product on the crystalline structure V, 7 S ^^(V). In terms of the mapping 
(f) it is defined as the pullback of the metric 5 on by the isomorphism i^^y 

(48) 7 = il,y9 , 

where ifp^y = d4){y) o ey and is the evaluation map ([T]) from V to TyM. The 
energy per unit mass 6(7) defines the thermodynamic stress vr, an element 
in {S2(V))', by 

where V = V{'j) is the volume per unit mass, related to the mass density fj, 
by 

1 

^(7) = —r^ , X = (l){y) . 

In the following, we assume the entropy a to be constant, what is called 
isentropic. 

Since we are in the relativistic framework, e includes the contribution 
of the rest mass-energy, so 

e = c'^ + e' 

(in conventional units), where e' is what we would call energy per unit mass 
in the non-relativistic framework. Note that the additive constant does 
not affect the definition (|49p of the thermodynamic stress vr, which may 



32 



MECHANICS OF CRYSTALLINE SOLIDS 



equally well be written in the form 

de' 1 

Also, in view of (jSOp below, we have 

E = Mc^ + E' , 

where 

M = [ dfi^ , E'= [ e'(7)d^^. 
Jn Jn 

The additive constant McP does not affect variations of E, which are the 
same as those of E' . For this reason we shall not distinguish in the static 
theory e', E' from e, E, and we shall denote the former by the latter. 

Remark II. 1. The static theory would be formally identical if formulated 
wholly within the non-relativistic framework. What distinguishes the two 
theories is that the non-relativistic theory is applicable only when the charac- 
teristic speeds i]i \29\l are negligible in comparison to c, while the relativistic 
theory holds for any values of the r]i less than c. 



2 Boundary Value Problem and Legendre-Hadamard Con- 
dition 



2.1 Euler-Lagrange Equations 

We choose a basis Ei,...,En of V such that uj{Ei, . . . , En) = 1, where 
Lo is the volume form on V. We denote by jab = i{Ea-,Eb) the matrix 
representing the configuration in this frame. We then have 



uj-f = ^(7)0; = det 7 • oj , 
and for the total energy of a domain in the material manifold 

(50) E = j e{^)dfi^ , 

n 

where di-Li_j is the mass element on M induced by to: 

, ^ , ( d 9 

(51) dfi^ — , . . . , - — 

\dyi y dyn 

Let {uj^ : A = l,...,n) be the basis for V* which is dual to the basis 
{Ea : a = 1, . . . ,n) for V. The oj^ are 1-forms on Af such that uj^{Eb) = 6^. 





D. CHRISTODOULOU AND I. KAELIN 



33 



Since {Ej\^{y) : ^4 = 1, . . . , n) is a basis for TyM, given any Y £ TyM, there 
are coefficients : A = 1, . . . ,n such that Y = Y'^EAiy)- Then: 

oj^iY) = (Y^EBiy)) = Y^iu;^iEB)){y) = Y^5^ = Y^ . 
Setting Y = we obtain 



Y^ = u^ 



d 

QyCl 



and we have e"i(y) = Y'^Ea G V. Thus 



(52) 



d 



^a{y)EA. 



Let us denote by dett<;(y) the determinant of the n-dimensional matrix with 
entries uj^{y). We conclude from (i5T]l . ([52]) that 



(53) = detw(y)d"y. 

The first variation of the energy (j50p is: 



(54) E 



d_ 
dX 



S(7 + A7) 



A=0 



A=0 



Ml) 



where 7 G S'^(V) is the variation of the configuration 7. By definition of 
the thermodynamic stress ()49p . we conclude, in view of (j53p . 

^ = - j ^7r^^7ABVdet7detu;(y) (Ty . 

Denoting by ruab the pullback hy (p to M of the Euclidean metric g on J^, 
we have 



(55) 
(56) 
(57) 



m 



rriab 



lAB 



9 9 

d]fdy^^'^ 



where the n-dimensional matrix with entries E'^[y) is the reciprocal of the 
n-dimensional matrix with entries io^{y). Thus det£'(y) = (det a;(y))~^, 
and by (j57|) we have 

det7 = det m (det ii^)^ , 
a/ det 7 det uj = V det m . 
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For the stresses S on M and T = (j)^S on M we have: 

(58) = -K^^E^E^^, 

(59) = S'^'^^. 

We will now formulate the Euler-Lagrange equations on both the material 
manifold and on space. We calculate for the variation of ^ab from (I57p 
using linear coordinates in E^ (i.e. dkgij = 0) 

The first variation of the energy thus reads 

/I r f)rr.j 

-S'^'rhabVd^ry = -J S'''g,, — ^Vd^d^y, 

Q n 
where is the variation of x*. Hence (1601) is 



(61) -/^"'^-■f^l^"'^-' 

n 

where d/i^ is the volume form on associated to the metric m. By the 
divergence theorem, (|61|) becomes 

S'''gi,-^Ma±' dAm + j x'VaQ'id^lm , 



(63) Qf = ^'^'5ra 



an n 
where Qi are the vectorfields with components 

dx^ 

m 

and V is the covariant derivative operator on associated to m. Thus, in 
local coordinates on A/", 



In ()62p . are the components of the outward unit normal to with 
respect to m, Ma = niabM^, and dAm is the area element of dfl associated 
to the metric induced by m. So M^ are the components of the covectorfield 
along d^l whose null space is the tangent plane to d^l at each point. 

The Euler-Lagrange equations are obtained by requiring E = for varia- 
tions x^ which are supported in and vanish on dO,. In view of (I62p . ()63p . 
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, the Euler-Lagrange equations are: 
d 



(65) — [Qydetinj =0 on . 



Note that with Qf defined in (j63p we have using (j56p 

(66) '^^W = ^-'^'h'^ = ^^^'^' = '^- 

Consider 



The first term on the right-hand side of (j67p is, by (|65p. (j66p . 



Thus (|65p is equivalent to: 

^3x* dy°-dy'' 



(68) V5,»=(^,— 



We shall show now that 

(69) "'^ "'^ 



dy"-dy^ ' 



so that the factor in parenthesis in (168p vanishes. In fact, (I69p is equivalent 
to: 

Qyd Q2^i 



dx^ dy°-dy^ 



which, since rucd^j = 9ij^^^ is in turn equivalent to 

(70) Tah.c = 9ij 



QyC QyaQyb 

m 

Now we have by the definition of the Christoffel symbols Tab,c and m from 
™ _ 1 / dnigc diTibc _ dmab \ _ dx^ d'^x' 

i a6,c - „ o,6 + r,„,a Q„,c - 9ij- 



2 \ (9y" dy^ J dy"-dy^ ' 

Thus (f70]) is indeed verified. We conclude that the Euler-Lagrange equations 
()65p are equivalent to: 



(71) VaSt = , or VaS'^" = on . 
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Remark II. 2. T*-' being the push-forward of S*"^ by (j) and rriab being the 
pull-back of gij by cf), the last equations (f7j[ j are in turn equivalent to 

(72) ViT^ = , or diT^ = on ^{n) 

in linear coordinates. 



2.2 Boundary Conditions 



Let now equations (j65p hold. Requiring again E = 0, this time for arbitrary 
variations x, we obtain from the first term of (j62p 

jab, 

an 



(73) 0= / S-'gi,^Ma±'dAr 



Here the covectorfield along di} with components Ma is defined by 

MaY"" = : Vy G Tj^aO , 

together with the condition MaY'' > for ah Y £ TyQ,, y G dfl, such that Y 
points to the exterior of and the normalization condition (m~^)"'' MaM^ = 



1. As (I73p is to hold for arbitrary variations i*, we obtain the boundary 
conditions 

(74) s''''g,~Ma = ^ ,or S^^Aff, = on 50 . 

Let X' = ^Y\ Setting M, = '^N^ we have: 

(75) NiX' = : G T.^^{dQ) . 

In fact, the covectorfield along (j){dO,) with components A^i is defined by (j75p . 
together with the condition that 

iViX* > : VA G r^.(/>(^^) , 

a; G (j){di}) = d{(f){Q)), such that X points to the exterior of and the 

normalization condition (ff"^)*"* NiNj = 1. In other words, := {g~^y^ Nj 
are the components of the outward unit normal to d{(f){0,)) in E^. In view of 
the fact that S"-^^ = T'^^, the boundary conditions ([74]) are equivalent 
to ^ 

(76) T'^Ni = on <P{dn) = d(p{n) , 

which means that no forces are acting on the boundary d(j){Q). 
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2.3 Legendre-Hadamard Condition 



The Legendre-Hadamard condition in the static case is just the second part 
of the hyperbohcity condition stated in (126p . see [4] for details. Here, it is 
formulated in terms of |^(y) = Va- I* requires for Ca,^b £ T*J\f,r]\r]^ G 
T^i^y-^M that 

1 9^6 

(77) ■^—-—j^aS.bvW ■ positive for ^, / . 

We differentiate e with respect to v^: 

de de d^AB 
dvi dlAB dvi 



where jab = ^A^Bdij'^a'^b^ hence 



-^f- - '^E{AEB)9ijv'a ■ 



{E^^E^^ denotes symmetrization in A and B) and thus 

We differentiate (fTSl) once more with respect to to get 
Therefore, 

1 3^ G 1 Qg cP" g 

C^^) 7 ^ - ^ A aihriW = \r]fiAiB + 7^ T. VCVA^B^D , 

4 dvl^dvl 2 djAB ojabOJcd 

where |ryp = Qijrfrj^ , £,a = -E'l'^a ^^"^ ^c* = E'^^^cdn'^f ■ Thus the Legendre- 
Hadamard condition reads: 

e 1 



At Tab the Legendre-Hadamard condition reduces to 



a2 



QiabQicd 



r]cr]AiBiD > (r/, ^ / 0) , 



because the first term on the right-hand side of (j79p (containing ^^^j—) van- 
ishes. Note that this condition is satisfied by virtue of the Postulate il. II that 
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o 



e has a strict minimum at 7, since this implies 



for any non-zero symmetric aA_B and, in particular, for 

CTAB = ^AVB + VA^B ■ 
In the general case, we define r]^ by rj^E'^v^^ = rf and we have 

\ri? = 9^,r1W = g^,vlviE\EU\'' = lABV^'v'' ■ 




Additionally, 





3 Uniform Distribution of Edge Dislocations in 2d 

The case of a uniform distribution of 2-dimensional edge dislocations is re- 
alized by the material manifold M being the affine group, see II. 1.11 The 

o 

corresponding metric n from (jlOp is isometric to the metric of the hyper- 
bolic plane. 

In general, we consider the total energy E ()50p of a domain Q in the 
material manifold M. In agreement with Postulate II.ll we assume that the 



energy per unit mass e has a strict minimum at 7- The symmetry of the 
problem motivates the choice of an isotropic energy function e that is a 

o 

symmetric function of the eigenvalues of m = (p*g relative to n. Therefore, 
in this model case, the crystalline structure is eliminated in favor of the 
Riemannian manifold {Af, n) . 

3.1 Isotropic Energy Function 

Concerning our toy energy function, we make the following basic choice: 



o 



(80) 




k=l 
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that satisfies e(Ai = 1, . . . , A„ = 1) = 0. Thus we have a strict minimum of 

the energy density at 7 =7- Note that, in ([SU]) . we are subtracting the rest 
mass energy, something which does not affect the variation of E, see also 
the discussion at the end of Section [H 

Therefore, 
(81) 



n 



<7) = ^X^(A|.-2Afe + l) = ^tr7'-tr7+^ = ^ ^13-^2^^^+^' 

k=l A,B=1 A=l 

Here and in the fohowing the basis {E^ ■ A = 1,... ,n) for V is chosen 

o o 

to be orthonormal relative to 7, hence 7a_b= ^ab- For such a choice to 
be compatible with the condition u}{Ei, . . . ,En) = 1, we must choose the 

o 

physical unit of mass so that the mass density corresponding to 7 is equal 
to 1. 

In case that the metric of the target space is Euclidean, gij = 5ij, we find 
from dSH) 

where the dot denotes the Euclidean inner product. For the total energy E 
(1501) of a domain 0, CC M, where uj coincides with the volume form on V 



(82) 7AB = E%E'—j-6,, = EIE_ 



induced by 7, we have d/x^j = d/io (the volume form on N corresponding to 

the Riemannian metric n), and we obtain for the total energy E expressed 
in terms of the components of 7, 

(83) ^= I \\iZ 71b-X;7aa + ^ I d//o. 

ii \ A,B=1 A=i J 

In the case n = 2, the invariants of a linear mapping are, respectively, 
the trace and the determinant. Expressed in terms of the eigenvalues Ai, A2, 
they read: 

T = Ai + A2 = tr m , 
6 = A1A2 = det m . 

We have: 

The toy energy ([8OI) thus reads 



e(Ai, A2) = I ((Ai - 1)2 + (Ai - 1)2) = 1 (r2 - 2r - 25 + 2) 
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3.2 Expansion of the Hyperbolic Metric 

We shall perform a dilation of the hyperbolic plane in order to be able to 
establish an existence result for a fixed domain, from which we can deduce 
an analogous result for energy minimizing mappings from a suitably small 
domain in the original hyperbolic plane to the Euclidean plane. 

We start with the usual metric, of constant curvature —1 for the hyper- 
bolic plane H, expressed in polar normal coordinates {R,9): 

ds]j = dR^ + sinh^ RdO'^ . 

Given a parameter e G M"'', we dilate H hy a. factor e~~^ obtaining Hf,. The 
corresponding metric, of curvature — e^, is: 

ds]j^ = [dR^ + sinh^ RdO"^) = dr^ + F^{r)de^ , 

where r = £~^R and 

Fe(r) = sinh (er) = (^er + ^ {erf + ^ {erf + . . . 

{r,9) are polar normal coordinates in ifg. The expression for F^{r) thus 
reads: 

(84) F|(r) = r'(^l + ^{erf + ^{erf + ...). 

The main point in the expansion (|84p is that F^{r) is a real analytic function 
converging to for e — t- 0, corresponding to the transition of the negatively 
curved hyperbolic plane to the flat Euclidean plane. 

In order to simplify the calculations, we go back to rectangular coordi- 
nates. We denote rectangular normal coordinates in if^ by (y* : z = 1, 2) 

= r cos , 
if = r sin 6 . 

Then 

dsl^ = [dy^f + {dy^f + e'{^- + ^^{erf + . . .) r'dO^ . 

Since 

d^ = -^ i-y^dy' + y'dy'') , 
denoting by ds'^ the Euclidean metric 

d4 = W + 
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the metric of takes the form: 

(85) = dsl + e'(^l + |(er)2 + . . .) {-y^dy' + yUy'f . 

The crucial point in the above expression for the metric h of is that 
it takes the form 

(86) hab = Sab + e^fab , 

where fab are analytic functions in {y^,y'^) and e^. In fact, 

(87) {-yW + y'dy^f = {r^5ab - y^y") dy'^dy' , 
hence, bv ([85]) . (|86]) and (|871) . 

(88) fab = Q + ^{erf + . . .) {rHab - y^y") • 
That is, 

fab = f{e\^) {r''6ab - y^y") , where fix) = 1 + Ax + . . . . 

4 Uniform Distribution of Screw Dislocations in 3d 

The case of a uniform distribution of (3-dimensional) screw dislocations is 
realized when the material manifold M is the Heisenberg group, see 11.1.21 

o 

The corresponding metric n from (jlip is isometric to the metric of a homo- 
geneous but anisotropic space. 

4.1 Anisotropic Energy Function 

The symmetry of the problem in this case motivates the choice of an anisotropic 
energy function e, since the dislocation lines have a preferred direction. 
Therefore, in this model case, the crystalline structure cannot be eliminated 

in favor of the Riemannian manifold {J\f, n) , in contrast to the 2-dimensional 
case previously discussed. 

The question of the right choice of energy per unit mass for the Heisenberg 
group is investigated in the following. Consider the orthogonal transforma- 
tion O of V, given by (X', Y' , Z') = 0{X, Y, Z) with 

X' = cosO ■ X + sme - Y , 

(89) Y' = -sine ■ X + cos9 -Y , 

z' = z. 
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The commutation relations [X, Y] = Z, [X, Z] = [Y, Z] = are preserved 
by this transformation since 

[X',Y'] = cos'^e[X,Y]-sm^e[Y,X] = {cos'^e + sm'^e)[X,Y] = Z = Z' , 
[X',Z'] = [Y',Z']=0. 

Consider the metric m = (j)*g. Since m is symmetric, we have generally 
six independent components in three dimensions. This also holds for 7, the 
corresponding inner product on the crystalline structure: 



7 



where the six components are 



01 P 



lAB = lAB A,B = l,2, 
Ba = -fA3 A =1,2, 
P = 733 • 

The energy density e depends on the variables (7, 0, /?). e must be invariant 
under orthogonal transformations ([89]) . that \se{^,9, p) = e{OjO, 06, p), for 
any O given by (j89p . |0p is such an invariant; together with the invariants 
tr7 and tr7^, which are the same as in the two-dimensional case, we have 
the following four invariants: 

(90) tij,tTf,\e\\p. 

Thus, starting with six variables, we have eliminated two and are left with 
an energy per unit mass of the form 

e = e{ni,p2, lOl"^ ,p) , 

O 

where /ii 2 are the eigenvalues of ^ab with respect to 7^^ in the Xy-plane 
iA,B = l,2). 

We may choose in particular: 

(91) e = i ((^1 - If + - I?) + f 1^1 V f (p - l)^ 

o 

which has a minimum at the identity jab =1 AB= ^AB- It is not surprising 
that the screw dislocations, which have a distinguished direction, break the 
isotropy and we are lead to an anisotropic energy of the form (I9ip . which is 
invariant under the transformation (1890 . 
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4.2 The Heisenberg Group as a Homogeneous Space 

Let 

^ dx ^ 5y ~'~ 5z ' dz 

be, as in Section [H the basis for the Lie algebra V of the Heisenberg group. 
We have the commutation relations 

(93) [X,Y] = Z , [X,Z] = [Y,Z] = 0. 
The basis of 1-forms (.^, ?]) dual to X, Y, Z is 

= dx , ( = dy , T] = dz — xdy . 

A left-invariant metric on the Heisenberg group is, up to isometries and an 
overall scale factor, given by: 

(94) n=^®^ + C®C + e^^V ® r] , 
that is: 

(95) ds"^ = dx'^ + dy"^ + e^f^ {dz - xdyf , 

where /3 is a real constant. The metric (|95|) is called a Bianchi type VH 
metric (for the classification into Bianchi types, see |12j ). It represents a 
homogeneous space which is not isotropic. 

o o 

The inner product 7 on V giving rise to the metric n is given in the basis 
{X,Y,Z) by: 

7(X,X) = 1{Y,Y) , l{Z,Z) = e^^, 

7 (X, y) = 7 (X, z) =7 (y, z) = . 

Thus defining 

(96) Ei=X , E2 = Y , E^ = e-^Z, 

o 

{Ea : a = 1, 2, 3) is an orthonormal basis for 7: 

(97) 1ab--=HEa,Eb) = 5ab. 
The dual basis is: 

(98) oj^ =i , ^2 = C , ^ e^77 . 

o 

In the following we shall denote the metric n by h. Thus, in terms of the 
basis {uj^ : A = 1, 2, 3) we have: 



h = y^^bj^ ® 
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Now (M, h) being a homogeneous space, any point in M may be taken 
as the origin. Consider a local coordinate system (y" : a = 1, 2, 3) on 
with origin a given point. Since in the Euclidean space we may set up a 
rectangular coordinate system (x* : i = 1,2, 3), a mapping id is then defined 
in the domain of the coordinate system on J\f by: 

id{y\y^y^) = {y\y^y^)€E^, 

that is, id is expressed in the respective coordinate systems as the identity 
mapping. Thus m, the pullback to M by id of the Euclidean metric 5ijdx^ (g) 
dx^ is 

m = dahdy"" ® dy'' 

in the local coordinates {y"" : a = 1,2,3) on M. The following problem 
then arises in connection with the analytic method to be presented below: 
given any point in M, find a local coordinate system {y'^ : a = 1,2,3) on 
with origin the given point such that jab, the components at y of the 
corresponding inner product on V in the basis {E^ : ^ = 1, 2, 3), that is: 

7AB{y) = mME\{y)E%{y) = 5abE'X{y)E%{y) 

satisfy: 

iAB{y)-iAB=0{\y?)- 

o 

(Recall that 1ab= ^ab-) The solution is simply to set up Riemannian nor- 
mal coordinates on [M, h) with origin the given point. For, the components 
of h in such coordinates satisfy: 

hab{y) = 5ab + 0{\y\'^) . 

Since 

hab{y)E\{y)E\,{y)=5AB 

identically, then 

lAB{y) - Sab = (Sab - hab{y)) E\{y)E],{y) = 0(|y|2) , 
as required. 

5 Scaling Properties 

5.1 The General Isotropic Case 

Let N = d M" and C E^ be the domain and the target space, 

respectively, of the mapping (p with the same properties as in ()47p . We fix 
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coordinates in M and also fix an origin in E"". We may assume that the 
coordinate origin in J\f is mapped by (p to the origin in E"". Consider 



ky) = (f ) , 



where / > is a given positive constant. If [/ is a domain in E"', we denote 
by lU the domain {Ix : x £ U}. Similarly for a domain in M". The domain 
of (/> is 17 = l^l, and 

The change from (p to (f) induces a change from m to m (and similarly for 

the inverses) as follows 

(99) 



1=1 



dya Qyb ""V/y ' ^ ! K^J X , yi 



For, since = Icj)"^ (j), we have 

V > Qya Qya / Qya 

We consider here the case of an isotropic energy function. We have seen 
that in this case the crystalline structure V can be eliminated in favor of a 
Riemannian metric h. 

Let us define a new metric /i on A/" by hab{y) = hah (j)- Note that h is 
not isometric to h. In the case of the toy energy function ()8ip . we have (see 
(fnill below), 



det h 

and similarly with /i, m, S replaced by /i, m, 5. Using (|99p we then obtain: 
(101) S''\y) = S" 



I. 

Now the equations ()7ip transform as follows. From the definition of the 
covariant derivative, we have 

(102) ViS"^^ = ^ + + ' 

m 

where F^^ are the Christoffel symbols with respect to m defined by 

™ 1 / f dnihfi dm..,] dmh^. 

(103) = -(,n-)- + - 
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as well as analogous expressions for V^S"^ and F^^. From (jM]) and (jlOip we 
then obtain 

(104) -^{y) = i^rA[l) and -rrriy) 



dy'i I dy<^ \IJ dy^ I dy^ \l 

and, consequently, from (jlOSp . 

(105) TUy) = ft (f 

We conclude, using (fTn2]) with (fTTM]) and (fT05]l . 

V,~S''\y) = \ (Vb5«^) (I 

So, if ij^ is a solution relative to h and </> is a solution relative to h and 
the tangent plane to Vl = 10. at the point ly G 90 being parallel to the 
tangent plane to dVt at y relative to the linear structure of M" D $7. 

If : A = 1, . . . , n) is an orthonormal frame field for h then [Ea '■ A = 
1, . . . , n), with Ea defined by 

E\{y) = E\ (I 

is an orthonormal frame field for h. Also, using a formula similar to ()105p 
with h in the role of m, we deduce 

i?,«,,(y) = l-^Rl, (I 

It follows that, if Kji is the sectional curvature of h corresponding to the 
plane 11 at y and Kf^ is the sectional curvature of h corresponding to the 
plane IT at y = Siy = /y, where II = dSiiJi) {Si the scaling map), then 

= r^Kn . 



5.2 The 2d Case 



Consider the two-dimensional case. The metric of the hyperbolic plane of 
curvature — is given in Riemannian normal coordinates by: 

Kb = Sab + ^'^fabiy) , 
fabiy) = f{e^r^)lab{y), 

lab{y) = Sabv^ - y^y" , 

and f{z) is an entire function with /(O) = |. 
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We now consider the metric h, the components of which in the above 
coordinates are: habiy) = hab{y/l)- We have: 

Kb{y) = 5ab + e^fab (y) , 

therefore, 

habiy) = 5ab + e^f (e^?-^) hbiy) = Sab + e^fab{y) , 

where e = e/l and fabiv) is fabiv) with e replaced by e. While the curvature 
K of h is — e^, the curvature K of h is 




Given a domain Oi, we shall show in Part III below that there is a ei > 
such that we can solve equations (j7ip for all < e < £i. Let (pi^s be 
the solution corresponding to ili and £ G (0,ei). ili is a domain in the 
hyperbolic plane of curvature — e^. We define (pi^e/i by 

(106) cPi^e/liy) = l^l,e (y) • 

This is the solution corresponding to the domain ^li = l^li in the hyperbolic 
plane of curvature —e^/t^. 

Choosing then / = e we have from ()106p 

^u(y) = ^'AM(y) =^0u(y) 

a solution of the problem for the domain Vti := IVL\ = eVii in the hyperbolic 
plane of curvature —e^/t^ = —1, that is, the standard hyperbolic plane. In 
conclusion, once we have a solution for the domain r^i and curvature — e^, 
we automatically have a solution for the smaller (rescaled) domain 0; = eVtx 
and curvature —1. 

Consider the stress T*-' and its rescaled version T*-' at the respective points 
in Euclidean space. From (jlOip and (jlOOp we conclude: 



Qya Qyb 

^ab,^.,d ^\y) d~^{y) ^ ry_\ d<P^{y/l) d<t>^{y/l) 

Qya Qyb \l) Qya Qyb 



_ rpij 
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The rescaled stress T at the rescaled point 4>{y) is thus the same as the 
original stress T at the point (j){y/l). 

5.3 The 3d Case 

Let {Ea : ^ = 1, . . . , 3) be the vectorfields (j96p that satisfy the commutation 
relations [Ei,E2] = e^E^, [Ei,Es] = [^2, ^3] = and (o;^ : ^ = 1, . . . , 3) 
the dual 1-forms (j98p . so that 

3 

(107) h=Y,^^^^^ > hAB=h{EA,EB)=5AB- 

A=l 

We introduce Riemannian normal coordinates (y" : a = 1,2,3) for /i at a 
given point in AA, which we take as the origin. The components of the metric 
h in these coordinates are of the form: 

(108) Kb{y) = Sab + eab{y) , eabiy) = 0{\yf) . 

Let id be the identity map defined in Section 14.21 and m the pullback by id 
of the Euclidean metric 5ijdx^ (8> dx^ . Then the components niab of m in the 
coordinates (y" : a = 1,2,3) are simply niab = ^ab and the corresponding 
inner product on V, which depends on y, is given by: 

lAB{y) ■■= j{y){EA,EB) 

= mabiy)E%{y)EUy) 
= 5abE\{y)El{y) 

= {Kb{y) - eabiy)) E%{y)EUy) 

= ^AB - eAB{y) , 

where eAB{y) = eab{y)E'X{y)E%{y) = 0{\y\^). 

We now dilate the Heisenberg group metric h by the factor / > 1, i.e. we 
set: 

h = l^h , h{EA,EB)=l^h{EA,EB) = f5AB. 
Define now Ea = 1~~^Ea, so 

h{EA,EB) = 5ab , 

that is (Ea : a = 1, 2, 3) is an orthonormal frame field relative to h. We 

denote by 7 the corresponding inner product on V, i.e. j{Ea,Eb) = Sab- 
Remark that the commutation relations of the frame field {Ea : A = 1,2,3) 
read 

[E,,E2]=r^E3 , [Ei,E3] = [E2,E3]=0. 
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Now: 

h = fh = fhab{y)dy'' ® dy^ . 
Set y"- = ly"-, then h = hab{y)dy°' ® y^ and hence 

(109) haM = Kb 

Consider the geodesic ray through the origin of the coordinate system y°': 

y'' = X't , y'' = \H where = IX" . 

For the Christoffel symbols of the metric h with respect to the coordinates 
y°' we have: 

r^,(At)A^A= = 0. 

But from (fT05]) . 

C(S) = irt (f) , 

hence: 

(110) f^,(At)A''A" = yr^,(Ai)/^A''A^ = . 

Thus (y" : a = 1,2,3) are Riemannian normal coordinates for h. 

We remark that the mapping id defined in Section 14.21 depends on the 
choice of local coordinates in A/". Let us denote by id the mapping associated 
to the coordinates {y°- : a = 1, 2, 3), reserving the notation id for the mapping 
associated to the original coordinates (y" : a = 1,2,3). Let fh be the 
pullback of the Euclidean metric 6ijdx^ dx^ by id. Then the components 
rhab of m in the coordinates (y° : a = 1, 2, 3) are again simply ifiab = 6ab- 

Now, by (fToHD and (fT09]) . 

hab{y) = Sab + eab{y) , 

where 

eab{y) = eab 

Therefore: 

iabiy) = O {l~^\y\^) , 
and we obtain, in analogy with the above: 

lAB{y) ■■= j{y){EA,EB) 

= mab{y)E\{y)&B{y) 

= SabEmEUy) 

= {hab{y) - ~eab{y)) E'X{y)E%{y) 
= Sab - ^Asiy) , 
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where eAB{y) = ^ab{y)E'\{y)E^{y) = 0{l~'^ \y\^)- Here E\ are the compo- 
nents of the vectorfield Ea in the coordinate system (y" : a = 1, 2, 3). 

Let now be a fixed domain in the y coordinates containing the origin. 
Let (/) be a solution of the boundary value problem 

(111) / Vfe5'^' =0 : in 0, 

\ S'^^Mb =0 : on aO, 

such that (j) takes the y coordinate origin in N to the x coordinate origin in 
E'^. In (fTTTD S"^ = ^'^^E'X&j^ and, by definition, 



/ det7 ,^g = _2- 



Vdet7 

Since h = hah{y)dy°' ® dy^ and = we have 

h = r^h = r^hab{y)dy'' (S) dy'' = hab{ly)dy'' ® dy^ , 

hence hab{y) = hab{ly)- We now consider the rescaled (smaller) domain 
^ = {y = y ^ ^} and define the mapping 

(^{y) = r^^{iy) -.yy^n. 

The pullback by (j) of the Euclidean metric is m, the components of which 
in the y"" coordinates are 

^^ab{y) = 5ij-^{y)^{y) . 

Since 



we have mabiy) = rhab{y)- 



Hence, the corresponding inner product on V, '^AB{y) = 7iy){EA, Eb) in 
terms of the original basis {Ea : j4 = 1, 2, 3), is: 

lABiy) = mab{y)E\{y)E'>B{y) = mab{y)E\{y)E%{y) = jAsiy) • 

For, 

Eliy)lL = l~'E%{y)-'' 



dy" ' 

hence E'^(y) = E'^(y). For the stress tensors in both coordinate systems we 
have vr^^(7(y)) = t^^^ {"yiy)) because intrinsically 6(7) = 6(7). Therefore: 

S^\y) = ^^^{^{y))EmE'B{y) = ^^''{l{y))E'X{y)E%{y) = S^\y) , 
which shows that the stress is scaling invariant, as expected. 
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Remark II. 3. The energy per unit mass e is a function of the configuration 
7, an inner product on V. However, when representing e as a function of the 
components 'jab of 'j in a basis {Ea : A = 1, . . . ,n), which is orthonormal 

o o 

relative to 1 , e in this representation depends indirectly on 7- 

Consider the equations in (jllip . Starting from riiabiy) = ^abiv), 
drhab,., ,^idrhab,.. .^^dniab , . 

the Christoffel symbols transform accordingly 

m m 

(112) vis) = i~'Kb{y)- 

Since we have 

(113) -w^y^ = ^ 'W^'^' 

we finally obtain, using (jll2p . (jll3p . 

) = ^^y^^ nMs'^y) + rlMS'^'iy) = (VbS'^'j (y) ■ 

Therefore, once we have a solution for the equation in the y coordinates, we 
immediately obtain the solution for the original equation and the boundary 
condition is also satisfied since Ma{y) = Ma{y). 
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Part III The Analysis of Equilibrium 

Configurations 

1 Stress Tensor 

We will now show that the stress tensor S on the material manifold M is 
given by 



(114) J^S^'' = -2-,^' 



det n 



dm„ 



Recall first the definition (j49p of the thermodynamic stress vr on the crys- 
talline structure V 



(115) J^^AB^_^de 



det 7 



djAB 



where we have made use of the definition of the volume ^^(7) = .1^^^ in 

y dct7 

o 

the frame {Ei, . . . , En) satisfying uj{Ei, . . . , En) = 1, 7yi_B= ^AB (so in fact 
det 7= 1). 

Now, from (|57p we have 

de de d-fAB de . 



dniab d-fAB drUab d-fAB ^ ' 

hence, with E{y) = uj{y)~^, 

de de A PI 

d-fAB dniab 

and finally, using ([58]) . (jll5p . (I116P and the equality 

, det m det 7 

(117) = ^, 

det n det 7 

we obtain: 



' detm ^ / det7 ^g ^ 9e ^ 9e ^ g 

det n V det 7 ^^^^ ^"^"^ 



a ^6 1 



which directly implies (jll4p . The equality (jll7p is a special case of the 
following proposition which applies to the isotropic case. 

o 

Proposition III.l. The eigenvalues Ai, . . . , 0/7 with respect to 7 coin- 

o 

cide with the eigenvalues of m with respect to n. 
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Proof. We define A{y) G £(V, V) by 

nA{y)YuY2) = ^{y){Y,,Y2) : yYi,Y2eV. 
Hence A{y)EB = A'^{y)EA, where 

Ai{y) = {V')^'^7Bc{y). 

Thus Ai, . . . , A„ are the eigenvalues of A. 

Similarly, we define B{y) G C{TyM,TyM) by 

n {B{y)Yi^y,Y2,y) = m{y){Yi^y,Y2,y) : yYi^y,Y2,y G TyM . 
and using B{y) = B^{y) ^ , we want to express A{y) in terms of 

B{y). 



Since the evaluation map Cy is an isomorphism from V to TyM for each 
y G A/", we have 

A{y) = ey^ o B{y) o ey . 
Hence, using Esiy) = E''^{y) ^ , 

y y 

Ai{y)EA = A{y)EB = ey\B{y)EB{y)) 
= E'B{y)Bt{y)oj^{y)EA, 

we conclude 

A{y) = E{y)B{y)co{y) = Lo{y)-' B{y)io{y) , E{y) = Lo{y)-^ , 

i.e. the linear mappings A{y) and B{y) are conjugate and therefore Ai, . . . , A„ 
are also the eigenvalues of B{y). □ 

Let us now calculate the stress tensor 5°^ corresponding to the toy energy 
(j80|) . Here, we denote n by h, the metric of the hyperbolic plane that is the 
metric of the material manifold for a uniform distribution of elementary edge 
dislocations. Recall that 

e = ^((Ai-l)2 + (A2-l)2) = ^(A? + Ai-2(Ai + A2) + 2) 

= ^tr/,(m^) - tr/,m + 1 

= \ {h^^y {h"^)'"^ mabmcd - [h"^)"^ niah + 1 ■ 
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Then: 

&°'--^aS;-^ ('.-)'V,.-(.-r) . 

or, 

(119) S'^>^ = 2^^{h~T{h~'f{K,-m^,) . 

The physical interpretation of (|119p is the following: 

i) If m is smaller than h then the stress is positive, 
ii) if m is larger than h then the stress is negative. 

Recall that a quadratic form q = h — m on TyN is said to be positive 
(negative) if, for all v S TyM , 

q{v,v) > 0(< 0) : Vt; / 0. 



2 Setup and Method in 2d 

As was shown in Part I, Section ll.l.H the material manifold J\f for the case 
of a uniform distribution of edge dislocations in two dimensions is given by 
the affine group, and a left-invariant metric on M gives M the structure of 
the hyperbolic plane of curvature — e^. 

To solve the problem in this case, we fix an origin in //^ and set up 
Riemannian normal coordinates (y" : a = 1,2) as in Section 14.21 of Part 
II. Let be any smooth bounded domain in these coordinates, containing 
the origin. Note that as e — >• 0, H^l^ tends to -f^olni where Hq = E is the 
Euclidean plane. We also choose an origin and set up rectangular coordinates 
(x* : i = 1, 2) in E. An identity mapping 

fl20) ^^'-^^ ^ ^ 

is then defined as in Section 14.21 of Part II. 

We now restrict the allowed mappings : C H/r E hy the following 
two requirements. First, 4> should map the origin in into the origin in E. 



Second, d(j){0) should map the vector ^ 



into a vector of the form A -A- 1 „ 

dx^ 10 



for some A > 0. By virtue of this restriction, the identity mapping ()120p 
(restricted to 0) is for e = the unique minimizer of our toy energy e from 
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Remark III.l. The restriction is needed to ensure uniqueness. Otherwise, 
composition on the left with a rigid motion of E gives another minimizer. 
The appropriate restriction in the three dimensional case will be stated below 
in Section \3.1l Analogously, it can be formulated in any number of space 
dimensions. Thus uniqueness is ensured in general. The argument of Part 
I, Section \3.^ applies with the hyperbolic plane Hf, and the Euclidean plane 
replaced by the n- dimensional hyperbolic space and the n- dimensional 
Euclidean space. This is as long as the toy energy Ii80\) is considered. 



From ()7ip and ()74|) the system of partial differential equations and the 
corresponding boundary conditions for the static problem of a uniform dis- 
tribution of edge dislocations in n = 2 dimensions is of the form 

where 



m 



(121) M 



We linearize the equations at the identity mapping, which is a solution for 
e = 0. First, we solve the linearized problem using the theorem of Lax- 
Milgram (as in [2]). An iteration will then show that there exists a solution 
to the nonlinear problem for sufficiently small e, and, therefore, by the 
scaling argument of Section [5] of Part II, that there is a mapping from a 
rescaled domain in the standard hyperbolic plane Hi = H oi curvature 
— 1 to the Euclidean plane, (p : Q C H ^ E, satisfying the conditions of the 
problem. 

We set 4> = id+ip, where ■0 is a small deviation from the identity mapping. 
We have: 

Fo [id] = . 

In a neighborhood of the identity is of the form 

F,[ct>]=Fe[id]+DidFe-iJ + N,[il^], 

where A'£[^/^] is to leading order quadratic in -0. We denote by the lin- 
earized operator DidF^. Then Ef, [(/>] = reads 

(122) Le-^ = -Fe[id]-Ne[>P]. 

We will solve this by an iteration starting at -00 = 0. In the first step of the 
iteration we have the linearized equations 

(123) • 01 = -Fe [id] . 



The term —Ef, [id\ in ()122p . (I123p can be interpreted as a source term. 
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A possible approach to the problem is to study the iteration 

(124) • V'n+i = -Fe [id] - Ne [^Pn] 

and show that ipn converges to a solution tp = (j) — id oi (jl22p , provided that 
we choose e appropriately small. However, we follow a different approach. 

Adding Lq ■ ip on both sides of (jl22p yields 

Lo • V = - (i. - io) • V' - F,[id] - N,[iP] . 
What we actually do is to consider the iteration 

(125) Lo • Vn+l = - {Le -Lo)-iJn-F, [id] - iV, [i^n] ■ 

Note that while both (jl24p . ()125p are linear in the next iterate ipn+i, in 
(jl24p the operator which refers to acts on V'n+i whereas in ()125p the 
operator which refers to Hq = E acts on ipn+i- 

In (jl25p . Le — Lq is a pair of linear operators, a second order operator in 
fi and a first order operator on dO,. The coefficients of these operators are 
of order e^, which we may write symbolically in the form: 

(126) Le-Lo^d{h-5) = e^df , 
since, from ([86]). h — 6 = e^/. 

The iteration (I125P starts also with -00 = 0. Then, setting n = in (jl25p . 
we have: 

(127) Lo-i^i = -Fs[id], 

thus, taking into account the fact that det/i = 1 + O(e^), (1119P implies 
S"^f'(id) = -2^2 + O(e^). It follows that is of order e^. This is 

step one of the iteration, the linear level. 

For n > 1, taking the difference between (|125p and the same with n 
replaced by n — 1, we obtain 

Lq ■ {lpn+l -Ipn) = - {Le - Lq) ■ {ipn " V'n-l) " {Ne[(pn] " Ne['lpn-l]) ■ 

In view of (jl26p and the fact that A'^j['i/'] is to leading order quadratic in 
V', for sufficiently small e, contraction will hold and ipn will converge to a 
solution -0 of (jl22p . Thus 4> = id + ip solves the problem (|12ip . 

Consider now the equations in ()12ip . To analyze the problem at the 
linearized level, we consider the variation iriab of the metric ruab at the 
identity mapping id. Setting 

(128) ^' = y' + , 
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and recalling that 



we obtain: 
(129) 

d 

rUab 



rriabiy) 



Qya Qyh 



ds 



51 +S 



dtp' 



dtp'' dtp" 

Qya Qyb 



Hence, niab is the Lie derivative of the metric 5ab on Hf, (the pullback by id 
of the metric 5ij on E) with respect to the vectorfield tp"'-^ on (which 

the push- forward by id takes to the vectorfield tp' on E) . 



3 Analogous Geometric Linear Problem 

In the work on the stability of the Minkowski space-time [6], a linear geo- 
metric problem was studied analogous to the linear problems in ()124p and 
(dlS]). 

Let {M,g) be a compact Riemannian manifold with boundary dM and 
X a vectorfield on M. We set 

(130) 7r = Cxg , 7r,j = V,Xj + VjX^ , X, = gijX\ 

the Lie derivative of the metric g along X, a symmetric 2-covariant tensor- 
field on M. The variation of tt with respect to X reads: 

(131) ^ij = ViXj + VjXi . 

We consider the problem of free minimization of the action integral 

(132) A = [ (h7T\l + p^X,) dfl, - I t'X, dflglg,, . 

Jai \^ J JdM 

Here p is a given vectorfield on M and r is a given vectorfield along dM. In 
mechanical terms p is the body force and r the boundary force. 

The first variation of (I132p . using (I13ip . is 

(133) A= I U^VjXi + p'X,)dpg- [ T^Xidpglg^. 

Jm ^ ' JdM 

For variations of Xi which vanish near the boundary we have 

Jm 

and requiring A = for such variations yields the Euler-Lagrange equations: 

(134) VjiT'^=p' :mM. 
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Requiring then ^ = for arbitrary variations: 

r'Xi dflglg^j = / {n'mj - r') X, dflglg^j = , 

JdM 

yields the boundary conditions: 

(135) n'^Nj = t' -.ondM. 

The equations (I134p . together with the boundary conditions (jlSSp . corre- 
spond to the linearized boundary value problem ()126p corresponding to a 
crystalline solid with a uniform distribution of dislocations in equilibrium if 
we identify {Vt,m,id* 5) with {M,TT,g), 6 being the Euclidean metric. 

We proceed in showing self adjointness of the above operators in the 
following sense. Let y be a vectorfield on M, a = Cyg, i.e. 

(136) =ViYj+VjYi = aji , Yi = gijYK 
Then we have by repeated partial integration using (jl30p . 

{Y, V • tt)l2(^m) ^^)l2(oa/) = V • a)^2(^M) N)^2(^qm) ■ 

Suppose now that y is a Killing field, i.e. a = Cyg = 0. Then we have by 
(fTMp . ([135]) and (fT36]) 

[ Yip'= [ y,V,7r*^ = -i / aij7T'^+ [ Yym,= [ Y^r', 

Jm Jm ^ Jm JdM JdM 

since aij = on M. Thus the integrability condition reads: 

(137) [ Y,p'= [ YiT\ 

Jm JdM 

This condition guarantees the existence of a solution for the boundary value 
problem (|134p , (|135p by the theorem of Lax-Milgram. 

3.1 Uniqueness of the Solution 

In fact, the solution is unique up to an additive Killing field. For, if we 
take two solutions Xi and X2 of ()134p with ()135p . their difference X = 
Xi — X2 satisfies the homogeneous equation of ()134p . i.e. /o = with zero 
boundary conditions, ()135p for r = 0. Therefore, we have, setting vr = Cxg, 
Vj {n'^Xi) = n'^VjX,, 

A = \[ \^\' = \l ^''^jX. = l[ V,(M,). 
^ Jm ^ Jm ^ Jm 

Using Gauss's theorem we obtain: 

A = \\ V,(^^^X.) = ^ / (7r^^iV,)X. = 0. 
^ Jm ^ JdM ^ — V — ' 

=0 
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It follows that vr = 0, and thus X = Xi — X2 is a Killing field, i.e. the 
solutions Xi,X2 only differ by a Killing field. 

For our problem, where (M, g) is isometric to a domain in the Euclidean 
plane, recalling the restriction 



d(/.(0) • 



m 

d 



0. 



_d_ 



we obtain, setting 
the conditions 



cP\x) = y' + sX\y) , 



sX\0) 

f)Vi 



0, 

X{s)6{ . 



Taking then the derivative with respect to s at s = yields the linearized 
conditions: 



(138) X\0) = 0, 

BX'' 

(139) ^(0) = Mi, 
where /x = A(0). The second of the above conditions is 

0. 



9yi 

Substituting the general form of a Killing field, 



X' 



a A 



-or- 



it follows from the first condition (|138p that /3* = 0, and from the second 
condition (|139p that a\ = 0, hence a*- = 0, i.e. X = 0. So the solution of 
the linearized problem is in fact unique. 

We remark that in three dimensions, one must add, for uniqueness, the 
condition that 

dm ■ i 

dy^ 



d 



and a^L. 

dx^ 10 



is a vector at the origin contained in the plane spanned by ^p- 1 ^ 
This can always be arranged by a suitable rotation in three-dimensional 
Euclidean space. Thus 



#(0) • 



_d_ 

9y2 



Ai 



_d_ 



+ A, 



d 
dx"^ 
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At the linearized level, this additional condition gives: 
(140) __(0) = ^i5i+^25^, 

where ni = Ai(0), H2 = A2(0). Setting i = 3 in (1140p . we have: 



while the conditions from ()139p read 

_(0) = ^(0)=0. 

Hence, if X is a Killing field, 

= af = and 02 = , 

or, taking into account that also /3* = (i = 1,...,3), we conclude that 
X = 0. 



4 The Linear Problem 

The unknown of the problem being the mapping (p : il. ^ E, the pullback 
metric m = (j)*5 depends on (p according to: 



Since mab{id) = 6ab, and (j) = id + ip, i.e. = + ip^{y), we have: 

Q^b Q^a 



rriab = Sab + mab + fJ-abiip] , where m« 



Qya Qyb 



and liabW\ is quadratic in -0. It follows from ()125p . together with ()126p . that 
ip = O(e^) at the linear level, whence ifP' = 0{e^). Furthermore, from the 
expansion of the hyperbolic metric in rectangular coordinates we have 

Kb = Sab + + C>(e^) , 

where lab = fab\e=0- 

Linearizing the operator Ff, from (I12ip at the identity mapping, LqiJ) = 
DidFo ■ ip, we find: 



Loip 



-2dbrnab ■ m Q , 
-2mabMh : on 50 - 
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with thab as above. Moreover, for = id we have, niab = ^ab, Vb = 
Hence: 

(141) ^S'^'iid) = Vd^ {h-')"' {h-'Y'' {ha - ma) = eXb + O(e') , 
where we have made use of 



Fe{id) 



Vdet/i = ^/det{5 + eH) = 1 + 0{e^) . 

Therefore, 

dbS'^^id) = 2e^dbl^^ + O(e^) : in n , 
S''^{id)Mb = 2eH''^Mb + O(e^) : on dVi . 

Thus, dropping terms of O(e^) in Fs{id) which come from the O(e^) terms in 
(jl4ip . the Hnearized problem (|127p reduces to the boundary value problem 

f _a_ . ^2, 
(142) <^ 



(143) 



[rhab - £ lab) = : in Q, , 



{mab - e^lab) Mb =0 : on (90 . 

We consider the following problem analogous to (|142p : 

Vj (tt^J -a^^) =0 : in M, 
(vr^i - o-*J) Nj =0 : on dM , 

where iTij = {Cxg)ij, cfij is a given symmetric 2-covariant tensorfield on M, 
A'j- is a covector whose null space is the tangent plane T^^M at x G dM, and 
A^* = {g~^y^ Nj is the corresponding outer unit normal vector to dM. 

To obtain the solution of (jl43p . we find a vectorfield X such that tt = Cj^g 
satisfies 

(144) - a'^) Nj = : on dM . 

We define X' = X — X and n' = Cx'g = Cxg — C-^g = n — n, to obtain 
from ()144p and the boundary conditions of ()143p : 

{Tr'fNj = {tt'^ - n'^) Nj = {n'^ - a'^) Nj = : on dM . 

Defining the vectorfield p in M by = Vj (cr*-^ — vr*-') , the problem then 
reduces to: 

Vj{TT'y^ = p' : in M, 

{ir'ymj =0 : on aM. 

This is of the same form as (jl43p . but with zero boundary conditions. To 
see whether this problem has a solution, we need to check the integrability 
condition (jl37p (i.e. orthogonality to the Killing fields in the sense). Here, 
the boundary terms vanish since we have zero boundary conditions and it 
remains for us to show that 

Cip'dpg = , 

M 
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for all Killing fields We have: 

Jm Jm 

= f ^^ [a'' - ^'') N, df,g\g^,j - [ V,C, {a'^ - 7^^^) d/i, = , 
JdM Jm 

where the first term is zero due to (|144p (that is, the choice of X), and the 
second one vanishes by virtue of the fact that ^ is a Killing field, 

V -e, {a'^ - TT^^) = i (Vi^, + Vj^i) {a'^ - n'^) = . 

By applying the Lax-Milgram theorem as in [2], we conclude that there is 
a solution X to the generalized linear problem which is unique up to an 
additive Killing field. Thus, also the linear case of the original problem 
(jl42p . viewed as a special case of the above, has a unique solution up to an 
additive Euclidean Killing field. 



5 The Nonlinear Case 



We will show that the nonlinear system (I12ip is solvable in Q under a certain 
smallness assumption on the parameter e. Let us first state the nonlinear 
case of the original problem (P), again. 



(145) (P) 



m 

V^S"^ =0 : m n, 
S^'^Mi, =0 : on d^, 



where 8°"^ is given by ()119p and M}, is a covector whose null space is the 
tangent plane TyQ. at y G dQ. The orientation may be defined by MfjX"" > 0, 
whenever is a vector pointing to the exterior of 0. M"" = (m~^)"^Mh 
is the outer normal to dil. The strategy for solving the problem (P) is 
the following. We set up an iteration scheme, where the first step is the 
linearized problem. By the result for problem (|143p . the linearized problem 
is solvable because the integrability condition is automatically satisfied. The 
integrability condition of the iteration can then be satisfied by applying a 
doping technique similar to the one in [9]. 



5.1 Iteration 



We first study the iteration scheme. For the analogous (generalized) problem 
(AP) we have from ()125p . ()145p (now written in terms of the coordinates y"" 
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on Q, and vr denoting the linearized metric m): 
(146) (^P)<^ -Pn ■ m 



K\i-<^)Mb =0 : on 50, 



where 



Now we set V'n+i = V'n+i ~ V'n+i, where the -i/'-part satisfies the boundary- 
conditions, i.e. 

(147) (<^+i-af)M, = 0, 

and we have a modified problem {AP') with zero boundary conditions 

, f <tl . • r. 

i TT^+iMb =0 : on af]. 

The integrability condition, which yields existence of solutions of the prob- 
lem by the Lax-Milgram theorem [7], now reads: 

That is: 

(148) / epl- I «'Mfe = 0, 

Jn Jan 

for every Killing field ^ of a background Euclidean metric on ft. This met- 
ric is id* 6, where 6 is the Euclidean metric of E. The y"" are rectangular 
coordinates on Q relative to this metric. In particular, is a unit covector 
relative to this Euclidean metric and the integral Jq^ is taken with respect 
to the measure on corresponding to the metric (arc length if dim il. = 2) 
induced on dil. by this Euclidean metric on Q. In the above, we have made 
use of (jl47p and the properties of ^ as a Killing field. In particular, 

,o,„,b I oa „a „b , „a I „b 



/-a „a„,b I oa „a „b . „a / „o „c 



and thus the contraction of the symmetric tensors tt^Xi, respectively a[ 
with TTTT vanishes. In concl 
step of the iteration is ()148p 



n ' 

with vanishes. In conclusion, the integrability condition at the n + 1 
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5.2 Killing Fields and Doping Technique 



In n dimensional Euclidean space we have the following linearly independent 
Killing fields: 

i) n translations, 
ii) "^"2 rotations. 

So the space of Killing fields on E"' is dimensional, where 

2 2 

Let {S^A ■ A = 1, . . . , N) be a basis for the space of Killing fields. In the spirit 
of [9] , we have to modify p, that is the inhomogeneity on the right-hand side 
of ()146p . This technique is called doping. We replace p by 

P = P + ^ caU , 

A 

and require in accordance with ()148p that 

(149) [^^.p'=[ ^j^.au : VA = l,...,iV, 

JQ Jan 

where a^j = a'^^Mi,. Since 

U<B = Mab 



is positive definite, we obtain a linear system of equations for the coefficients 
CA ■ ^ = 1, . . . , as follows: 



u-p'= U-P+ uly2cB^B]= u 

Jn Jn J Jan 



■<TM ■■ VA=l,...,iV^ 



and thus 

MabCb = / iA- CTM - I iA- P =■■ CFA ■■ yA=l,...,N . 

^ Jan Jn 

The system YIb ^'^ab^b = (^A can be solved, and there is a solution 

(150) CA = Y.{M-^)j,j,aB, 

B 

Mab being positive definite, hence non-singular. 
We reformulate the problem {AP) as follows: 

(151) (AP)l ■ 

I {<''+i-<')Mb =0 : on dn. 
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Writing V'^+i = X^ = /J 
d I dX^ 

(152) 



/ a -.ab 



afMb = t", (fTKTD is of the form 



I ax^ 

dyb \ dy^ "I" dy'' 



dx'' I 

^ dy" 



dy°- 



P 



in Q. , 
on 50 . 



Under the restriction discussed above, which if X is a Euchdean Kihing field 
forces X to vanish identically, the linear system ()152p has no kernel. The 
following estimate then holds (see [l]) 

(153) ||X||^^^,(^)<c{||p'||^^(o) + ||T||H,^^/^(e^)} , c = c{n). 

Applying this estimate to ()15ip and taking e suitably small we can prove con- 
traction of the sequence (V'n) in Hs+2{^) for s > §, Hs{^) and Hg_i/2{dO,) 
being under this condition Hilbert algebras. 



Cn,A 

in Hs{n) pl 
in /f,+i/2(9f^) at 
in i?s+i(f^) tt: ■ 
in Hs{^) 



OA 
Ca 

pa 



ab 
n+l 

rn 



ab 
ab 



> for n 



oo , 



and, in the limit n — )• oo we deduce in terms of the modified inhomogeneity 
P' 

(154) 



a- 



ab 



dyb 



P 




in il. , 
on Of] . 



We have: 
(155) 



o^' + Y.caCa, 



and CA is given in terms of a a by (jl50p . From (|145p . (jl54p and (|155p . we 
find 



(156) 



'W 



We have to show that 
Proposition III. 2. 

(157) X'' = Y,CAeA = ^ 



1,. 



in , 
on dil . 



if e is suitably small. 
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Proof. Let be a Killing field of 5ab- We have: 

(158) / 5abex'd^y= [ mabCX'd^fim + 0{e^\X\J, 

Jn Jn 

where we say that a real-valued function Q{e,X) is 0(e^ if there is 

a constant C such that 

\Q{E,X)\<Ce'\X\^ . 

(jl58p holds because ruab = 6ab + O(e^). Consider the vectorfield C on J], the 
push-forward of which by (p to the Euclidean plane E coincides with ^ 



Qya 



Then is a Killing field of the metric m = (j)*6 on Q. Hence: 
(159) 



VbCaS''''dfim+ [ CaS'^^'Mb d^lm\^n > 
JdQ. 



where Qa = mabC' ■ The integral on Q. in (|159p vanishes because 



ab 



The integral on 90 vanishes by virtue of the boundary condition in (|156p . 
Then, since 

^'^ = + 0(e2) ^ 

we deduce: 



(160) 
However 



/ SabCX^'d^y = [ mabCX'dfim + 0{e^ \X\J = 0{e^ \X\ 



a = 1, . . . ,n . 



and since Killing fields are analytic functions, we have on a bounded domain 

n-. 



(161) 

From (fTBOl) . 
(162) 



\X\ 



C max Ic^l . 
A 



< Ce^ \X\ 



U-Xd^y 
But from the definition of X in (jl57p , 

/ ^A-Xd''y = Y^ / UiBCB d^y = ^ Mabcb = (J A , 
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and hence, from (|162|) . 

(163) \aA\<Ce^\X\^. 
On the other hand, we see from (jl50p that 

(164) maxlc^l <Cmax|(T^| , 

A A 

where C = ||M-^||. Finahy, using (fT6T]l . (fTBl) and (fT63]) . we obtain 

l^loo = ^1 k^l ^ ^2 max |cj^| < Cae^ l-'^loo > 

which imphes X = for e sufficiently small (Cse^ < 1). This finishes the 
proof. □ 

Remark III. 2. Concerning the regularity of the solution for C°° domain 
(smooth boundary) the solution is also C°°. For, ^ G Hs+2, s > n/2 implies 
ip G //s+3(r2). Therefore, by induction, ip G Hk{ft) for every k. 
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